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1.0 Introduction 


The research activities supported by this contract include the evaluation 
of CCD detectors for x-ray astronomy applications, contributions to the development 
of an Imaging gas scintillation proportional counter, the evaluation of certain 
metal oxide crystals as potential radiation detectors, optical observations and 
searches for x-ray sources discovered by the HEAO-1 A2 experiment, and theoretical 
modeling of non equilibrium ionization structure of supernova remnants. 

2.0 Detector Development 

2.1 Evaluation of CCD detectors for x-ray astronomy 

A program of CCD evaluation for x-ray astronomy was initiated under a President' 
Fund at CalTech (PF173) and the Jet Propulsion Laboratory in Pasadena, California. 

The purpose of the Investigation was to irradiate existing CCD devices being used 
in optical wavelengths with x-rays to determine their performance characteristics 
as x-ray detectors. The funding was quite limited and only allowed for the testing 
of a few existing devices. The initial tests were not very encouraging, but did 
show that the optical CCD detectors were sensitive to x-rays. The energy resolution 
was very poor and the efficiency was low. In order to continue this work, we 
proposed to use somo grant money NAGW-108 as well as requesting an additional 
President's Fund for one more year of effort. The extra support was still not 
adequate to produce CCD's more optimally designed for x-ray detectors, but did 
permit more testing and some electronic modifications in the readout processing. 

One immediate benefit co both the Space Telescope Project and the Galileo 
Project was that the x-ray data provided a good evaluation signal (which were 
hard tc obtain using visible light) to determine parameters for the optical chips. 
This situation occured because the x-ray interaction generated a precise amount of 
charge in the CCD which could be used to evaluate charge transfer characteristics 
in the silicon material of the chip. 



Improvements in signal processing of the charge transfer in the three phase 
clocked CCD f s resulted in a much improved performance. These results were pre- 
sented at GSFC on October 6, 1981 (Riegler et al 1981). A new virtual phase 
device was evaluated for this report as well. All CCD's were produced by Texas 
Instruments since they were the ones being used for Space telescope and the 
Galileo Project. 

Further testing has continued following the above report. Considerable 
improvements have been made as reported by Stern et al 1983. At this time the 
energy resolution and detection efficiency are close to the theoretical values for 
silicon. It is likely, that with the production of x-ray optimized CCD's which 
are now being fabricated at Texas Instruments, the resolution which is set by the 
energy required to produce an ion-electron pair and the Fano factor for silicon 
will be realized. 

A proposal was submitted on the basis of the above results and has been 
funded by NASA for making X-ray optimized devices and to fly one of these detectors 
on a sounding rocket to evaluate background levels and detector performance eX 
the focus of a Wolter Type I x-ray telescope. 



3-0 Gas Scintillation Proportional Counter Observation of M87 

An Imaging Gas Scintillation Proportional Counter (IGSPC) was developed by 
the Los Alamos National Laboratory's Astrophysics branch under the direction of 
Dr. Doyle Evans. Dr. Richard Blake together with Dr. David Anderson undertook 
the detector development at LANL. A 15" diameter Wolter Type I mirror assembly 
(see Tuohy et al 1979) was combined with the IGSPC to form an imaging instrument 
with superior energy resolution and high detection efficiency. The energy 
resolution achieved for the flight instrument was 7.9% at 5.9 keV, with a position 
resolution of 1mm at this energy. Pulse rise time information was available for 
background rejection. The following sections describe the preliminary results from 
the rocket flight which took place on 21 June 1982. 

3.1 Scientific Rationale and Choice of Target* 

The IPC (Imaging Proportional Counter) flown on the Einstein Observatory 
prosessed good detection efficiency over the 0.2-40 keV energy band, moderate 
spatial resolution ^1 mm, poor energy resolution, 120% at 1.5 keV, and some 
background rejection efficiency ^50% of the cosmic rays. In order to improve on 
this capability for future x-ray astronomy missions, better position resolution, 
energy resolution and background rejection present themselves as possible areas for 
investigation. The gas scintillation detector had been developed over the previous 
decade, to where the possibility for a substantial gain in energy resolution could 
be made. The positioning of an event, however, posed a more difficult problem. 
Optical systems could be used to locate the light pulse in the scintillation chamber 
either as in an Anger Camera used with crystal scintillators, or by imaging of the 
scintillation light on an optical imaging detector. 

A clever solution to this problem was discovered by David Anderson, at 



LANL who showed chat by adding a small amount of TMAE to a proportional counter 
gas that it could be rendered sensitive to the UV light produced by scintillation 
in Xenon gas. This discovery made it possible to combine the superior energy 
resolution of the gas scintillation counter with the position resolution of the 
imaging proportional counter. Since this technique preserved the inherent good 
timing characteristics of the proportional counter, the rise time of the light 
pulse could be measured at the same time as the position and energy thereby 
providing the possibility of improved background rejection. 

The Einstein Observatory has obtained extensive observations of the brighter 
supernova remnants and clusters of galaxies, both classes of targets which are 
the only ones available for detailed observations from sounding rockets. In 
order to obtain improved Information over the Einstein results, it is necessary 
to make observations with either improved spectral and/or spatial characteristics. 
The IGSPC offers such a possibility in that the IGSPC has the spatial resolution 
of the IPC and the energy resolution of the SSS, thereby offering in one 
observation the capability of both instruments on the Einstein Observatory. The 
two brightest targets with adequate spatial extent and appropriately emitting 
spectral properties matching the combined telescope-detector 

response are Pup A and M87. The latter source is a little low in intensity, 
but not so low as to be Beverly photon statistics limited. 

3.2 The Instrument Performance and Initial Results 

The IGSPC was flown with the 15" Wolter Type I telescope from WSMR on 
21 June 1^82 using a boasted Black Brant (27.068 AH). The flight was completely 
successful, in that M37 was acquired for the entire flight of some 250 seconds 
above 130 km altitude. The aspect camera showed that M87 was held to witnin 
about one resolution element during the entire exposure, and positioned very 


near the center of the detector. 



A calibration source illuminated the detector with AtK x-rays both at the 

a 

beginning and end of the t light. The pulse height distribution is shown in 
Figure 2-1, indicating an energy resolution of about- 22%, a factor of 1.5 
better than a typical proportional counter resolution such as obtained on HEAO-1. 

The rise time distribution of events produced by the source provides 

a measure of the x-ray induced time spread. The cosmic ray and gamma ray Induced 
events come from the region outside of the telescope field of view. Their 
distribution of rise times is much broader, especially toward longer rise time. 

The rejection efficiency for this particular counter is about 70% at 1.5 keV. 

The x-ray rise time distribution was not found to be energy dependent between 
1.5 and 5.9 keV. 

The position resolution of the detector has been measured in the laboratory 

as a function of energy. At 5.9 keV the resolution is about 1mm. At AtK the 

oi 

resolution degrades to 2mm, becoming progressively worse as the energy is lowered. 
At this level of position resolution, the angular resolution of the experiment is 
dominated by the detector and not the mirrors which have about a 2 arc minute 
resolution (note that 1mm corresponds to 1.7 arc minutes in the focal plane). 

The rocket flight only recorded about 600 x-rcy events from M87. This is 
less than anticipated which was mostly caused by a high threshold on the iow 
energy discriminator to reject noise in the emplifiers used on the cathode wires 
for position sensing. 

About one third of the events originate from the cluster core where possible 
heavy element sedimentation and rapid cooling of the hot gas may be occuring. 
Detailed modeling and fitting of the data has not been completed at this time. 
Given the lov; statistical quality of the data, it is unlikely that no more than 
upper limits will be set to evolution in the cluster core region. 



Figure Captions 
Figure 1 

In flight calibration of the IGSPC prior to and following the M87 observation. 
The energy resolution is about 22%, FW HM. 

Figure 2 

The rise time distribution of events compared to the cosmic ray background. 

Figure 3 


The spectrum of the core, inner 4 arc minutes, and the halo of M87. Spectral 
modelling of this data has not been completed. 
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4.0 Non-Equilibrium X-ray Emission from Young Supernova Remnants 

Young Supernova Remnants (SNR) emit strongly in the x-ray band of the 
spectrum. Numerous lines from a variety of elements have been observed in their 
emission spectrum (Winkler et al 1981, Becker et al 1980), Because tb^M; ends > 
lines are generated in a hot dynamic medium which is undergoing rapid expansion, 
the emission lines are produced in a non equilibrium condition (Moore and Ganrire 
1976, Itoh, 1978). Since young SNR's still contain a substantial fraction of 
the ejected material of the parent star relative to swept up interstellar gas 
(ISM), it is of great scientific interest to attempt to determine the actual 
abundance of elements in the SNR with the goal of determining something about the 
final stages of stellar evolution and the enrichment of the interstellar medium 
via supernova explosions. Detailed modeling is essential, since the plasma is 
far from collisional equilibrium. The following calculations are only a first 
step in this investigation since the detailed nature of the initial conditions 
for any given SNR such as the density distribution of the local ISM and the 
initial blast energies are never known with precision and certainly vary from 
remnant to remnant. For a detailed discussion of the hydrodynamics and cal- 
culational methods refer to the CalTech Thesis by the same title written by 
John Nugent Jr. 

4.1 Outline of the Model 

The hydrodynamics are calculated assuming a spherical supernova explosion 
of energy, E q , which expands into a uniform ISM with particle density, n^. The 
energy is assumed initially to be almost completely contained in the ballistic 
motions of the ejecta. The flow of the gas, except at shock fronts, is assumed 
to be adiabatic, so heat conduction and radiative cooling are neglected. An 
explicit set of finite-element equations representing the Lagrangian equations 
were used in the calculation. These equations, which were taken from Richtmyer 
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and Morton (1967^, included the artificial viscosity term, which acts as a 
dissipative mechanism at shock fronts. The number of spherical shells used for 
the ejecta was 13 and the number for the shocked ISM varied from 12 to 24. The 
number of shells was constrained by CPU memory and time considerations. Dis- 
continuities in density at the ejecta/ISM interface and at the jump at the main 
forward shock were smoothed slightly because of the finite number of shells 
resulting in some minor underestimate of the emission from these regions relative 
to the rest of the remnant. 

Thus, given E q , n^, and the mass of the ejected material, M fi , the mean tem- 
perature, T^, and the electron density, n e , can be calculated as a function of 
time and radius. The calculation of the electron temperature, T^, depends on 
the question of non-Coulomb processes at the shock front. Two extreme cases 
will be examined. The first is to assume that the electrons come across the 
shock front "cold" and are subsequently heated by Coulomb-like collisions with 
the ions. The second is to assume that some mechanism at the shock front y 

brings the electrons and ions into thermal equilibrium. Under those conditions 

T ■ T « T . 
e i m 

Once the history of T g and n g in a given shell is known, the fractional 
population of ionic state with charge, ?. of an element with nuclear Z, c*n be 
found by integrating 


n 

e 



*Z,z-l 


(i_ +r_ )X_ +r_ . . X_ , 

Z,z Z,z Z,z Z,z+x Z,z+1 


where X_ is the fractional population, n i is the rate for an atom with 
tL , z e 6 y z 

nuclear Z in ionic state z to ionize to state z + 1 via electron collisions, and 

n r_ is the rate for an atom with nuclear charge Z and ionic state z to recom- 
e Z,z 

bine with an electron to form ionic state z-1. The fractional concentration is 


defined such that 
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The set of values {X^ X^ z+1^ 8 ' ia ^ hereafter be referred to as the 

ionization structure of the element with nuclear charge Z. 

Some pre- ionization must occur in front of the shock for a collisionless 
shock front in relatively low stata? of ionization. Hydrogen and helium are 
assumed to be initially stripped. The rest of the elements considered and their 
given initial state are: C +4 , N +5 , Ne +5 , Mg +5 , Si +5 , S* 6 , Ar +6 , Ca* 6 , Fe* 6 , and 
Ni . Because of the relatively large ionization rates that still occur at these 
early stages of ionization, the transient effects caused by this ad hoc choice of 
initial conditions should be eliminated in % 10-20 years. 

The final step is to use the electron temperature, electron density, and ioni- 
zation structure to calculate the spectrum from each element in each shell. A 
modified version of the Raymond and Smith plasma code (1977) is used. The 
collisional excitation rates for He-like ions were changed to those calculated by 
Pradhan, Norcross, and Hummer (1981), which, unlike earlier calculations included 
effects of autoion izing resonances. Also, a number of processes that produce Kot 
nhotons from iso-sec j^nces with lower ionic charge than He-like ions have been 
added. The original Raymond and Smith code only considers innershell excitation 
for Li-like, Be-like, and occasionally B-like ions. This is sufficient in a CIE 
plasma because if the temperature is high enough to excite Ka emission from 
lower iso-sequences (i.e. if % E^) then the iso-sequences below B-like are in 
very low concentrations. This is not necessarily the case when the plasma is in 
non-ionization equilibrium (NIE) especially for the higher Z elements. For the 
densities and time scales applicable in the SNR problem, we felt it was adequate 
to augment the Ka processes only for Ar, Ca, Fe, and Ni. Innershell excitation 
to the L-shell has been added, then, for a.i iso-sequences between Li-like and 
F-like ions of these elements. In addition, for these elements, innershell 



ionization for all iso-sequences below Li-like and the process of producing K 
shell vacancies via di-electronic recombination has been added for iso-sequences 
between Li-like and O-like ions. Both these processes result in Ka emission. 
Excitation rates, branching ratios, and resulting photon energies for these 
innershell processes are taken from Gronenschild (1979 and references within 
(e.g. Mewe, Schrijver, and Sylwester (1980)). 

Parameterization of the model 

The input parameters to the model and the units we will always use are the 
initial blast energy, E q (units of lO^ 1 ergs), che density of the interstellar 

_3 

medium around the progenitor star, n Q (units of 1 cm ), the amount of ejected 

mass from the progenitor star, M g (units of Mg) , the age of the remnant, t^ 

(units of 1000 years), the abundance, relative to hydrogen, of elements in the 

ejecta (normalized to the measured solar value), the abundance, relative to 

hydrogen, of elements in the ISM (normalized to the measured solar value), and 

the distance to the remnant, D (units of kpc). We also will define a "filling 

factor", f, that gives the fraction of the solid angle of the remnant that is 

filled with shocked plasma. Finally, we also fit for the amount of neutral 

22 2 

hydrogen, N„ (units of 10 cm ) intervening between the remnant and the earth. 

H. 

This latter quantity parameterizes the energy dependent attenuation along the line 
of sight due to photoelectric absorption and scattering. The Fireman model (1974) 
was used for this calculation. 

Since the calculation of the spectrum for a choice of E^, u 9 M g , and the 
a^e, takes a relatively large amount of computer time, it is necessary, in 
order to do spectral fitting, to interpolate from a grid. By exploiting a 
scaling law, we need only to store the results on a three dimensional grid. As 
discussed by a number of authors a given set of parameters E^, n^, M^, t^, and D 
does not produce a unique spectrum (e.g. Mansfield and Salpeter, 1977). An 
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infinite set of those five parameters can be found that produce the same spectrum 
for the same abundance set. This degeneracy can be discussed in terms of a scaling 
parameter defined as 


a 



1/3 


as well as a number of "reduced" parameters: 


2,1/3 -2_ 

n - (En ) - a E - an, 

o o o o 


0 * * n -1 a" 2 M , 

E* 6 


and 


n \ 1/3 


ft 


- «- l t 3 . 


A fourth reduced parameter is the reduced emission volume, C which is defined as 

ev 

/n 2 dV 

c = - e -- 9 - 

ev 4 ttD 2 


and can be expressed as 


aT 'f 

* . ev - 

b- A 

ev 4-D 2 


where T ' is a function of n,6, and t and is defined as 
ev 


T ev' (n,a,T) H a_1 J * n e 2<iV * 
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13 -5 

Wc will quote C #v in units of 10 cm . 

Those four reduced parameter along the and the abundance sets of the 

ejects and ISM form the parameter a;>ace searched when fitting spectra. Once a 
point in this reduced space is determined to be the best fit, the next problem 
is to Invert from the fitted parameters (n.S.t, and C #v > to the physical super- 
nova parameters (E ,n ,M ,t,,f, and D). So far we have shown four equations for 
six unknowns. A fifth equation can be added that relates the angular radius of 

the shock front, 6, a measured quantity, to R ' , which is defined by 

s 

R ’(tf.t) ^ «~ l R 

S 3 

where R is the true sire of the shock radius and R ' is a function of t> and t. 
s s 

Since 0D - R * we can write 
s 

R ’ - a -1 0D . 
s 

This still leaves the svstem of equations unde rdete trained by one variable. Our 
approach will be to assume f to be unity and show how the deduced parameters 
scale with f. Taking that approach* the equation for C^ v and the equation for 
can be combined to yield the following expression for o: 

T 1 (n t s* • i )0 2 
ev , 

* • f 

AnR (t? * t )C 

s ev 

With the input parameters are given by 

2 2 -1 
E • M ■ ■ ot n “an 

o e 3 o 

and the distance is given by 

R * 



Since a is directly proportional to f, each parameter scales with f as it scales 
with a. 

We digress at this point to give some physical significance to this n-8-x 
space. As can be seen from its definition above, 6 can be thought of as an 
“ejected mass-like" axis and x can be thought of as an "age-like" axis, n can have 
a few different inte re tat ions* but for our discussion it is best described as a 
"density-like" axis. The 8x plane is where the shape of the pertinent hydro- 
dynamic variables , electron temperature and density as function of radius are 
determined. If we let X be the ratio of the radius to the shock radius, then the 
temperature profile, T(X) and the density profile divided by n* n *n(X), are 
functions of B and i only. Finally, a is measure of the mean x-ray surface 
brightness. 

Various regions in the B-x plane can be related *o the different phases of 
the SNR evolution. An Important quantity for discussing this demarcation is x 

s 

which is given by 

t - .16 iS 5/6 . 

s 

When T v x s the remnant is in the so-called "free expansion phase". During this 
phase the ejecta dominates the hydrodynamics of the remnant which expands at a 
constant rate which depending upon the initial velocity of the ejecta. At 
i ■ t , the mass of the shocked ISM begins to exceed the mass of the ejecta, and 
the hydrodynamic profiles begin to approach those of the Sedov similarity solu- 
tion for a point explosion (Sedov, 1959). The x-ray emission from the ejecta is 
signit leant for x v Beyond that time the contribution of the ejecta emission 

drops below a few percent of the total emission. 

Increasing B increases the emission volume of the ejecta component relative 
to the emission volume of the shocked ISM and decreases the ratio of temperature 



in the shocked ejecta to the temperature in the shocked ISM. For x < x g the 
temperature of the shock region is inversely proportional to B as well. 

Increasing t for a fixed 8 means a larger fraction of the x-ray emission comes 
from the shocked ISM as shock wave sweeps up more interstellar matter. The 
temperature of the shocked ISM drops monotonically with increasing x - slowly 
during the free expansion phase and like x during the Sedov phase. The 
temperature of the ejecta initially rises with x during the formation of the 
reverse shock but later, at the times of interest for this work, it also decreases 
with increasing x. For a fixed 8 and x,n is a scaling factor for the density 
profile, t\t is proportional to / n^dt which is the effective time for heating 
electrons in the "cold" electron model and the effective amount of time for 
evolving the ionization structure. Increasing n increases T^/T^ (for Coulomb 
heating model only) and results in the ionization structure of the heavy elements 
to be closer to CIE. 

Coir; irison with previous work 

Early efforts for modeling non-equilibrium SNR, x-ray spectra were made by 
Itoh (1977, 1979). He constructed models with and without electron-ion thermal 
equilibration behind the shock front, but n and 8 were both fixed in his papers. 
The lodel in Itoh’s 1979 paper was done with no ejecta at all. Concurrently with 
ti: . model developed for this thesis. Shull (1982) and Hamilton, Sarazin, and 
Chevalier (HSC; 1982) have also devised similar models. Shull (1982) calculated 
models only with n ® 1, with the no ejecta, and only with electron-ion thermal 
equilibrrtion assumed. HSC calculate models varying in n and x, for both 
elactron-ion post shock equilibration and for only Coulomb electron heating, but 
again no ejecta contribution was considered. This work differs from these other 
efforts in that it includes the reverse-shocked ejecta as well as the shocked ISM. 
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The models by Itoh (1979), HSC, and Shull rely exclusively on a Sedov similarity 
solution to provide the hydrodynamic substrate. This latter approach is equivalent 
to this work when t >> t • In this "far Sedov" limit the spectrum is completely 
independent of 8, and, as is done in Shull (1 (, '82) and HSC, t can be converted to 
a shock temperature as follows: 

T - 5.2t _6/;> keV . 
s 

We also note that the definition of n is slightly different in HSC. For comparison 

n jjn 

This thesis also is among the first efforts at fitting spectral data with NIE models 

where n varys. (Shull and HSC currently have work in progress.) 

These three other NIE models developed plus the one described here are nearly 

independent* At least two different methods for obtaining the hydrodynamic 

quantities, three or four different codes containing the atomic processes, and 

different codes for integrating the electron temperature and ionization structure 

were employed in these computer models. Despite these differences and many other 

idiosyncratic differences, all models seem to agree with each other where it is 

4 

possible to check. For software systems of this size (^ 10 lines of code) and 
subtle complexities (numerous) this is not a trivial statement. 



The Far Sedov limit 

The Sedov limit will be an important limit in this thesis and a convenient 
limit for discussing some of the properties of the process of mapping fitted 
parameters to physical parameters. In this limit, 




and 
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R ' a: 1.5 X 10 19 t 2/5 cm 
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From these relations, a for a Sedov limit model is given as 




a - .If 


arcmin 


'13 


where ® . is the angular radius in units minutes of arc and C,, is C in units 

arcmin 13 ev 

of 10^ cm f, again, is a filling factor. The distance is given by 


D » 1 B 


OT 


2/5 


e 


kpc . 


arcmin 


As can be seen from above, in the Sedov limit the final set of parameters 
scale quite differently with n,t, and f: 


e --nV'V t,* n V' 5 f D.„ 2 ,‘ /5 f n . ,fV 2 / V l 

O J o 


The higher the powers to which n*t, and f are raised the more severe the error 
propagation. So E q is the least well determined while n o is the best determined. 
This leads us to Nugent's Corollary of Murphy’s Law of SNR modeling: "The more 

interesting a parameter, the more poorly it is constrained". While the blast 
energy, an important parameter for such issues as stellar evolution and 
nucleosynthesis of heavy elements, is typically uncertain by an order of magnitude, 
we can, with dazzling precision, measure the amount of garbage around the 
progenitor star. 

In the way of a digression, the reader may be bothered that D is an output 
of the model as opposed to an input. In previous x-ray papers discussing SNR's, 



the parameters known were T (from which T could be inferred under the 

e s 

assumption that * T^) , the angular radius, and the x-ray flux. There were 
three equations relating these measured quantities to E^, n^, t , and D. Since 
there were four "unknowns", one had to be fixed so that the other three could be 
determined. Unless the age was known from historical records, the distance was 
usually the best determined "unknown". Heuristlcally speaking, with the advent 
of the NIE models another quantity, namely / n e dt, is measured from the amount 
of ionization heating and, possibly, from the amount of Coulomb heating. This 
gives a fourth equation so that all parameters can be determined. 

Most distance estimates to SNR's rely on the empirical I - D technique using 
radio measurements. In some cases, additional estimates have come from analyzing 
historical records of the peak apparent magnitude of the supernovae event, or the 
progenitor star was a member of a certain star cluster. The use of NIE models in 
fitting x-ray spectral data represents an independent measure of the distance 
that arises from a non-empirical, physical model. 

Caveats 

Before proceeding to discuss the fits of the NIE model to x-ray spectral 
data, a few caveats are in order. First, despite the touting above, some care 
must be taken in accepting the inferred values of E q , n^, t, M^, and D. Their 
rather fragile derivation could be destroyed by a number of systematic problems. 
Some of the problems that are suspected or are known to exist in SNR's include: 
inhomogenities may exist in the ISM; the ejecta/ISM interface is Rayleigh-Taylor 
unstable; and shocks might not form over the entire solid angle of the remnant. 

Heat conduction may be another important effect, especially if dense ejecta or 
ISM clouds are heated to shock temperatures. All of these may have a profound 
effect on the emergent spectrum. In addition to these hydrodynamic considerations. 



systematic uncertainties in the atomic physics as well as detector modeling could 
also swamp any of the statistical uncertainties. 

In order to put the NIE model in perspective, a brief review of the history of 
x-ray spectral modeling of SNR emission is helpful. At first, multi-temperature 
thermal bremsstrahlung models were fit to the data with the caveat that line 
emission and non-equilibrium effects were being ignored. Then as data and 
theoretical work Improved multi-temperature CIE models that included the line 
emission were employed. The lack of line emission was no longer a problem but 
non-equilibrium conditions along with the number of problems outlined above still 
were. These models were taken to be parametrizations of the problem, but it was 
always the expectation of the observer that the parameters derived were indicative 
of the real remnant. With the introduction of the NIE model, the lack of con- 
sideration of the non-equilibrium effects has been removed. In light of the 
caveats given above, it is safest to still view the NIE model as a parametrization. 
To a rough approximation the model is representing a plasma with a characteristic 
electron temperature and an effective ionization time (possibly a second plasma 
if the ejecta is a significant component to the emission). As was discussed above, 
before NIE models the angular size, x-ray flux, and implied temperature were used 
to derive remnant parameters. The NIE model has added another measurement, the 
effective ionization time, and otherwise improved the derivation procedure; 
however, the basic idea is the same as before. 



4.2 The Observation of MSH-14-63 (RCW86) 

rlSH 14-63 is the most likely candidate for the supernova recorded by the 
Chinese in A.D. 185 (Clark and Stephenson, 1977). While MSH 14-63 is the name 
of the radio source, an alias for the remnant is RCW 86, which is the designation 
of the optical filaments. Following the discovery of MSH 14-63 (Hill, 1967), other 
early radio maps clearly resolved a limb brightened shell with an angular radius 
of ^20 arcmin (Caswell et al., 1975). Westerlund (1969a) suggested that the rem- 
nant was associated with an OB association. If this were the case, a Type II 
supernova would be indicated on the basis of stellar type. By using spectropho- 
tometric studies, Westerlund determined the distance to the association to be 
2.5 kpc. This measurement is in agreement with the measurement of the distance 
made by the I - D technique. 

The first x-ray observation of MSH 14-63 was made by Naranan et al. (1977) 

in the range 0.5 to 2.5 keV using an experiment on the Apollo-Soyuz mission. The 

observation was taken in the range 0.5 to 2.5 keV. The spectrum in the range of 

1 to 30 keV obtained from the 0S0-7 spacecraft was reported by Winkler (1978). 

Winkler found that the combined spectrum of the two observations could not be 

explained by a simple thermal spectrum with a single temperature. Rather, two 

components each at different temperatures were needed to adequately explain the 

data. The temperature of the high-temperature component, which is presumed to 

be produced by the shocked ISM, was measured to be greater than 5 keV. The low 

temperature component, measured to be 0.22 keV, was attributed to the supernova 

ejecta which had been heated by inward propagating shock (McKee, 1974; Gull, 1973, 

1975). Using this reverse-shock model, Winkler determined an ejected mass of 

greater than 5 M- is required. Winkler concluded that the x-ray spectral data 
w 

supported the idea that MSH 14-63 was formed from a Type II supernova located in 
the OB association studied by Westerlund. 



Because of Increased spectral resolution and sensitivity of observations, 
most recent attempts to explain the SNR spectral data have employed more 
sophisticated models to predict the x-ray spectrum from hot, optically thin plas- 
mas than the simple exponential spectrum used by Winkler. These models explicitly 
calculate the contributions from each element in the plasma to the continuum 
emission as well as the the line emission (Tucker and Korea; 1971, Kato, 1976; 
Raymond and Smith, 1977; Shull, 1981). For a large range of conditions of the 
plasma line emission from heave elements (Z 6) dominates the spectrum, and as 
was discussed in the last chapter, a simplifying assumption that is often made 
is that the plasma is in collisional ionization equilibrium (CIE) . 

In this chapter we will present results of both CIE and NIE models with data 
collected from MSH 14-63 by the A-2 experiment on the HEAO-1 spacecraft. 

Observations 

The observations were made by two different detectors of the A-2 experiment 
(Rothschild et al., 1979) on board the HEAO-1 satellite. The low energy pulse 
height data (.5-3 keV) were obtained using the on-axis LED detector while the 
spacecraft was in a scanning mode. The high energy pulse height data (2-15 keV) 
were taken during a pointed observations using the MED detector. 

Each detector had two separate co-aligned fields of view. Roughly, the 
dimensions of these two fields were 1?5 x 3° full width at half maximum (FWHM) 
and 3° x 3° (FWHM)* The data from both of these fields were combined to obtain 
the final data set used for the analysis. The dimensions of all the fields of 
view of the detector were larger than the dimensions of MSH 14-63 (diameter of 
40 arcmin) , so the data represent the integral over the entire surface brightness 


of the remnant. 



The LED detectors had a sensitivity to photons as low in energy as 0.18 keV; 
however, because of contamination from a Centauri (Nugent and Garmire, 1978) the 
lowest energy data were not used. Since interstellar absorption of the spectrum 
from MSH 14-63 is so great at energies below 0.5 keV, this was not a significant 
loss of information, a Centauri is measured to have a steeply falling spectrum 
(kT^ x 0.15 keV) and in broad band measurements in the energy range above 0.5 keV 
c* Centauri has a count rate £ .04 that of MSH 14-63 (Nugent et al.» 1982). These 
facts, combined with fact that the collimator transmission to at Centauri during 
the MSH 14-63 observations was % 0.2, gives us confidence that there has been no 
significant contamination of the pulse height data above 0.5 keV. 


CIE Fits 

We will first discuss fits of CIE models to the MSH 14-63 data-both a single 
temperature plasma and two-temperature plasma. There are three reasons for this 
digression: 

1) Because bremsstrahlung and radiative recombination from hydrogen and helium 
dominate the spectrum at photons energies £ 4.5 keV, the electron temperature 
and reduced emission volume alone determine the shape and magnitude of the 
spectrum. The determination of the temperature* as we will discuss, will 
eliminate the Coulomb heating models. 

2) Since only a two-termperature model was used by Winkler, we would like to 
apply the same model for purposes of comparison. 

3) We would like to demonstrate how poorly a single temperature CIE model fits 
the data, before proceeding with NIE model fits. 

Two-temperature CIE models assume two separate coronal plasmas of different 
temperatures. In addition to the temperatures and the reduced emission volumes 


of the two different components, other parameters are the elemental abundances 
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of the two components, and the interstellar absorption column, N„. 

H 

2 

The results of the fit are shown In Figure la. A reduced X of 1.6 for 60 

degrees of freedom Is obtained. The only elemental abundance that was allowed 

to vary was iron. The results of the fit are shown in Table 1. The parameter of 

particular Interest to this discussion is the higher temperature, .>.1 keV, which 

will be denoted, T , . , . As mentioned above, at photon energies £4.5 keV the 
6 » nxgn 

spectrum is dominated by continuum emission from hydrogen * helium and, hence, 

T u.r u is insensitive to NIE effects. 
e,hign 

An important result from this value of T , . , is that we exclude the whole 

e,high 

branch of NIE models which assume the electrons are heated by Coulomb collisions, 
Cox and Anderson (1982) , using analytic approximations to the Sedov solutions for 
the temperature and densities in the interior of the remnant and analytic approxi- 
mations for Coulomb heating, showed that the electron temperature is independent 

3/14 -4/7 

of radius for times shorter than ^3000E e , n . This ”• iateau** temperature 

ji o 

as a function of time is given by 



In terms of our fitting parameters this equation becomes 


T % 1.6 


-2/25 


n 2/5 keV 


As an example, Figure 2 shows the results of numerically modeling coulomb heating 
for n * 0.2 and t ■ 1.2. The two top graphs shou ^he density and mean temperature 
shown with an analytic Sedov solution. The third graph shows the electron 
temperature assuming Coulomb heating only. The dotted line which displays Cox 
and Anderson 1 s approximate solution demonstrates excellent agreement. 



Because is insensitive to the choice of parameters, an attempt to find a 

choice for D, E . and n that allows T * T . . and satisfies the constraints for 
o o e e,uigh 

angular radius and reduced emission volume produces an absurd result. To wit: 

D ^ 41 Mpc, E q ■ 3 x 10^ ergs, and n Q » 6 x 10 * cm \ We will safely ignore the 
possibility, then, that the electrons are "cold" after having cresses the shock 
front. 

Saying that there exists at the shock front some mechansim other than Coulomb 
interactions to couple the electrons to the ions does not prove taac the mechanism 
w?U be sufficient to bring the electrons and ions into thermal equilibrium. If, 
after passing through the shock, the electrons have obtained a significant 
fraction but not the total amount of their equilibrium energy, then the time scale 
for them to reach equilibrium would still be as above. The only difference would 
be that the electron temperatures would be higher than those derived by Cox and 
Anderson (1982). 

Before exploring the fits using the NIE models, we wish to emphasize that a 
single high temperature component does net fit the data at low photon energy ' 

Figure lb shows only the 5.1 keV component with a reduced emission volume as 
fitted in the tv^-temperature model. The flux at 1 keV differs from the predicted 
flux by a factor ^ 10. 

Review of Parameters 

Before discussing results of the model, we briefly review for the reader the 
definitions of the parameters used. The input parameters to the model and the 
units we will always use •'re the initial blast energy, E q (units of 10^ ergs), the 

-3 

density of the interstellar medium around the progenitor star, n Q (units of l cm ), 
the amount of ejected mass from the progenitor star, (units of Mg), the age of 
the remnant, (units of 1000 years), the abundance, relative to hydrogen, of 
elements in the ejecta (normalized to the measured solar value), the abundance, 
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relative to hydrogen, of elements in the ISM (normalized to the measured solar 
value), and the distance to the remnant, D (units of kpc). These parameters are 
mapped to more convenient set of parameters for purposes of the doing spectral 
fitting. There is a scaling parameter defined as 


a 



1/3 


and a amber of "reduced" parameters: 


2 . 1/3 - 2 _ 

n - (E n ) - a E - an , 

o o o o 


M . « 

6 - -rA “ n *a M , 
E e 

o 


and 


T 



1/3 


C 3 


-1 

a 


t 


3 ' 


A fourth reduced parameter is the reduced emission volume, C fiv , which is defined as 


C 

ev 


jn ^dV 

Sr 


n.i.iJ, and C are derived from fits to the spectral data, a is then deter- 
ev 

mined from the average x-iay surface brightness by 

1 c 

ev 

where V* is determined from n, and * action of solid angle of the 

remnant filled with shocked plr. ... 1 of f will nominally be assumed to 

equal unity. With a, the input parameters are given by 
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2 2 -1 
E * a n M ■ a n6 t- * ax n » a n 
o e Jo 

and the distance is given by 

R '(6,t)qi 

“--nr 

where R 1 is a constant determined from the fitted value of 6 and x* 
s 

When x ^ 16g^^ the remnant is in a Sedov phase and the shock temperature, 

T , can be given as a function of x by 
s 

T - 5.2 I " 6 ' 5 keV . 
s 

Non-Ionization Equilibrium Models/T^ ® 

We are going to assume for the balance of this paper that whatever mechanism 
is responsible for coupliug electrons to ions at the shock front will bring them 
to equilibrium on an ignorably short time scale. Under this assumption then we 
will always assume that T * T. - T . 

The parameter 6, as it turned out, is not a well constrained parameter for 
tais data set. Results of fitting the data with models that have relatively high 
values of 3 will be discussed below; however, we will assume, for simplicity and 
definitiveness, during most of the discussion of the MSH 14-63 data that 6*1. 
Since the values of that fit the data under this assumption are ^1.1, a deli- 
berate consequence of this choice of 3 is that we will be working in the Sedov 
limit. There are three advantages to this: 

1) Working in a regime of parameter space where emission from the shocked ISM 
dominates Che spectrum is deemed preferable to working in regime where the 
reversed shocked ejecta was dominant since the modeling of the latter is 


more controversial. 



2) Since the spectrum is Insensitive to B in this limit, the results of the fit 
are also Insensitive to the precise choice of 6, other than it is sufficiently 
low, 

3) Models that are in the Sedov limit make possible direct comparisons with the 
work of Shull (1982) and Hamilton et al. (1982). 

In addition to fitting for C fiv , N^, n* end t, we also allowed the abundances 

for Fe, S, and Si to vary. The abundances for the rest of the elements were fixed 

at solar values as given by Allen (1976), The results of the best fit are displayed 

in Table 2 along with the input parameters inferred from assuming the radius of 

the remnant to be 20 arcmin. Figure 3 shows a display of the best fit and an insert 

2 2 

showing the dependence of K on n and t. This X grid is calculated by fixing the 

values of n and t and letting the balance of the parameters vary. The inner con- 

2 

tour is for change of X from minimum of 4.6 which corresponds to the minimum 

volume 90% confidence region for 2 parameters (Lampton et al., 1976), For com- 

2 

parison, a contour for AX ■ 10 is also shown. The results of fits that are done 
at two extremes of the 90% region are shown in Table 2. 

Ue will now discuss various aspects of this fit: 

Elemental Abundances 

The over-abundance of the S and Fe was needed primarily to fit the lines 
at 2.4 and 6.7 keV. Because L shell emission from iron is a major contributor to 
the spectrum at ^1 keV, the reader may be concerned that the overabundance was 
required to match the flux in the LED data. To demonstrate that this is not the 
case and that a solar abundance composition can fit the data with the exception 
of the lines, we have also performed a fit where all the elemental abundances in 
the N1E model were fixed at solar, and lines were inserted, ad hoc, ,>t the Ka 
energies of S, Fe, and Ar (2.4, 3,0, and 6,7 keV, respectively) to compensate 
for the missinf flux at those photon energies and in effect neutralize the in- 
fluence of the lines in the fit. The results of the fit are given in Table 3 and 
the fit shown in Figure 4. 



The effect of changing the Fe abundance does have an e iect on other parameters. 
In this case* increasing n has the effect of compensating for the lowered eler eu:al 
abundances since the population of Fe XVII and Fe XVIII* both Important 1 keV 
L-shell emitters* Increase dramatically as n is raised from 0.13 to 0.3. 

Density* age* and 

The fit to the combined LED /MED data gives good agreement to other known 
properties of MSB 14-63. The density of 0.1 -0.3 is reasonable for the height of 
the remnant off the galactic plane (% 80 (D/2kpc)pc). The measured value of 
would imply an average hydrogen density along the line of sight to the remnant 
of '0.7 (2 kpc/D). This is also a reasonable value considering the angle of the 
remnant off the galactic plane. 

Most important lv, the fitted values of the age are in reasonable agreement 

with the age of 1800 years known from Chinese astronomical records. Alternatively* 

if one assumes an age of 1800 years* then the average x-ray surface brightness is 

consistent the measured value. This not a trivial statement since the average 
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x-rav surface brightness scales as n . a large fraction of the volume of the 
remnant cannot be filled with a plasma which had a much higher n than those 
fitted above. 

Fe K line 
a 

The Fe is an ;t tractive feature for study because it is prominent and 
measurement of the surrounding continuum is easy. The line is composed of a 
Mend of lines around ' 6.6 keV from many different ionic species of iron. There 
are two anomalies concerning the line in the NIL model: 

1) An over-abundance of iron must be assumed to explain the equivalent width of 
the feature. 

2 ) The centroid of the line feature in the data is apparently shifted to higher 
energies with respect to the predicted line feature. These shifts in the cen- 
troid are ' l‘>0 eV to lower energies. Doppler shifts can probably hi ruled 



out since the shock velocity implied from the x-ray temperature could at most 
produce a shift of ^AO eV if all the iron-line emitting material wa on the 
far side of the remnant. 

In a 5 keV plasma which is in CIE, the most prominent lines arise from H- 
like and, more so, He-like ions, but when the plasma is out of CIE the populatiow 
of these ions can be severly diminished. Given the ionization time, in 

MSH l*-63, predicted by the above NIE fits the prominent iso-sequences are S-like 
to F-like ions. The emission, which results following innershell ionizations of 
these ions, has photon energies of ^200 eV lower then the characteristic Ka radi- 
ation from He-like and H-like ions. At greater ionization times the ionization 
structure becomes more advanced. The characteristic energy of the line photons 
increases with increasing ionic charge. The efficiency of ions below C-like for 
producing is generally lower than the strong He-like ion by a factor of ^3. 

The C, B, and Be-like ions, on the other hand, produce comparable amounts of Ka 
from innershell excitations via dieltctronic recombination. 

To reconcile the observed iron feature in the MSH 14-63 data two requirements 

in some region are required. First, in order to produce Fe Ka emission at all, 

T £ 5 keV. Second, the effective ionization time would have to be £ 5 times 
e 

that of the models used in the fits abo^e. The latter requirement is necessary to 
start populating the ionic states at or above the C-like isosequence. The charac- 
teristic radiation from these ions is not only at higher photon energies, but 
again, these ions also have higher efficiencies for producing Ka. 

In order to quantify this discussion, we will perform NIE model fits on just 

the data ^4.5 keV. Using these data will avoid many assumptions involved with 

fitting the data at lower energies which are dominated by line emission. We will 

consider only four parameters: n> t, C , and the abundance of iron, A^. The 

amount of N , which only marginally affects the data at these energies was fixed 
H 



at the fitted value from above and the abundances of all other elements were fixed 


at solar. In way of review before proceeding, t is primarily responsible for the 
shape of the continuum, ht is proportional to /n g dt , is proportional to the 
equivalent width (strength) of the line, and C gv acts as an overall scale factor. 
t and C gv are well determined by the continuum and will not change by more than 
10% for this discussion. 

Figure 5a shows a fit done with n fixed at 0.2 and the balance of the parame- 
ters left to vary. The non-iron contributions to the spectra have been subtracted. 

As noted above, the model shows a distinct shift of 'v 150 eV to lower energies and 
an abundance enhancement of factor of 3.2 over solar is needed to explain the 
intensity of the feature. If we now unpin the value of n, a best fit is obtained 

with a value of n = 3.1. This fit is displayed in Figure 5b. Not only does the 

2 

line fit better at this value of n (The improvement in X was 'v 11), but it does 

so with an abundance much closer to solar. The 95% confidence region (2o) for n 

ranged from 1 to 6. (Even at n * 6 there are still significant deviations for CIE.) 

Finally the same region was fit with a CIE model. Now a single electron 

temperature, T, has replaced t. The fit displayed in Figure 5c, also shows a well 

centered line with near solar abundance; however, the CIE model predicts too much 

emission at ^ 8 keV. This line feature, which has been referred to as K8 (e.g. 

Pravdo and Smith, 1979), contains n®3 to n=l transitions from H-like and He-like 

Fe and Kot emission from Ni. The Fe component tends to dominate because Ni is less 

abundant than Fe by a factor of ''-IS. The ratio of K8 to Ka has been used by 

Pravdo and Smith as a diagnostic for CIE since it is very sensitive to the population 

2 

of the He and H-like ions in Fe. The X of the fit to data points between % 5.5 keV 
and ^ 9 keV increased by 12.1 from best fit NIE to best fit CIE. This increase, 
which is evidence against the plasma being in CIE, came primarily from contributions 


around the KS feature. 



Taking the best fit value of n from the Fe-line fits what can one say about 

2 

the remnant? First, because the surface brightness (or a) scales as n , we can 
rule out a model where plasma of this nature fills the entire volume of the rem- 
nant. For filling factors near unity, the distance would need to be 200 kpc and 
the age, 3 x 10^ years. Filling factors of ^1% would be adequate, although we 
know from x-ray images probably 50-75% of the remnant is filled with plasma. 

These x-ray images (Pisarski, private communication) show a limb brightened shell 
that has large variations in surface brightness as a function of angle around the 
ring - the brightest region behind in the SW. In addition there are noticeable 
deviations from a circular ring. There is some correlation for brighter regions 
to have a smaller radius of curvature. All of this suggests that there probably 
are density inhomogen i ties on various different size scales present in the remnant. 

Since the models calculated above assume a homogeneous ISM, they are what 
we would call "single-n" model. What the spectral data and the x-ray image seem 
to require is, in a sense, a "multi-n" model. Somewhat analogous to "two- 
temperature*' models, we would argue for low n's to explain the shape and flux of 
the low energy emission and high n's to explain the shape and flux of the high 
energy lines. 

lhe high-n components could for example be explained as dense cloudlets that 
are heated by conduction to temperatures of the post-shock region. Another 
possible soxrce of material in the remnant with high ionization time is the 
ejecta. The reverse shock cannot be the source of the Koi emission since the 
temperatures are not high enough (McKee, 1974; Gull, 1974, 1975). There have 
been two mechanisms suggested for mixing the ejecta into the hotter, snocked ISM 
regions - neither was appropriate for inclusion in one-dimensional hydrodynamic 
model. If the ejecta was initially in a smooth shell, the contact surface be- 
tween the shell and the shocked ISM would become Rayleijh-Taylor unstable 
(e.g. Gull 1975). Fingers of this dense ejecta material could then break away 



from the rest of the ejected material and be heated to temperatures comparable to 
the shocked ISM by heat conduction (Chevalier, 1975). Alternatively, soae of the 
ejecta could initially be in clumps to begin with (Chevalier 1977), I*, edition 

to being heated by shock waves due to pressure differences between the clump 
interior and the shocked ISM, the higher temperatures could, again, be obtained 
by heat conduction. 

Distance 

We turn our attention next to the question of the distance to MSB 14-63. As 
mentioned in the introduction, Westerlund (1969a) has measured the distance to an 
OB star association in the direction of MSH 14-63 by spectrophotometric technique. 
He finds a distance modulus of 12.0 + 0.2 which corresponds to a distance 2.5 kpc 
with an uncertainty of roughly 0.2 kpc. This is inconsistent with the distances 
related to the x-ray Measurements discussed here which argue for the distance 
being much closer (^1.6 kpc). 

Given that the age is known, the best method for discussing the limits on 
the distance by x-ray measurements is, again, to work only with data above 4.5 keV. 
Here the only assumptions are that the remnant can be approximated by a Sedov 
blast wave, that the electron temperature is equal to the mean temperature, that 
the age is 1800 years, and that the angular radius of the shock front is 20 arcmin. 
The distance is then given as 



The minimum volume 95% conficence region (2o) in t extends from 0.9 to 1.2. The 
lower limit of 0.9 (an upper limit on the shock temperature of ^6 keV) gives an 
upper limit on the distance of 1.7 kpc. For the x-ray determined distance to 
agr oe with Westerlund result would require t ^ 0.55 (plasma temperatures of 
12-15 keV) , which can be ruled out. It is, of course, possible that the progenitor 



for MSH 14*63 was not a member of Westerlund's OB association; however, we do 
point out that any deviations from the assumptions stated above would tend to 
make us underestimate the distance. 

First, even though the longest angular dimension of the x-ray emission is 40 
arcmin, the brightest regions * those dominating the x-ray spectra - appear to 
have a smaller radius of curvature, maybe as small as IS arcmin. Second, if, 
despite some non-Boulomb electron heating at the shock front, T e is still less 
than T^, than the shock temperature would be underestimated by our measurement 
of T x * Finally, the introduction of density inhomogenities into the remnant can 
alter the interpretation. If, for example, after expanding for a significant 
part of its lifetime, the shock in the SW part of MSH 14-63 moved into a denser 
region, the shock temperature would be lower than the radius would imply. 

Finally, while a temperature of 12-15 keV can be ruled out if the x-rays were all 
arising from a post-Sedov shock, a thin, hot component related to the postshock 
regions could exist that was being overwhelmed by cooler emission from denser 
regions. The observations lack the sensitivity to detect such a hot component 
under these conditions. 

A lower limit can be set for the distance by assuming that the remnant is 
still in the early free expansion phase. A lower limit of the kinetic temperature 
of the shock of 5 keV would imply a constant expansion velocity of 2100 km s \ 
Assuming, again, an age of 1800 years and an angular radius of 20 arcmin yields 
a lower limit on the distance of 700 pc. 

Emission from the ejecta 

The value of d has consequences for demonstrating the existence of the 

reverse-shock ejecta component and measuring the mass of the ejected material. 

As was mentioned above, 6 was a well constrained parameter. We fit the data 

with values of d ranging from 1 to 10 with no more than a change of roughly unity 
2 

in the value of X* - . At low values of d the ejecta was only a few percent of the 



total emission at ^ 1 keV and, thus, did not significantly effect the spectrum. 

At higher values of 8 the fraction of emission from the ejecta became dominant 
but the temperature of the ejecta decreases. The ejecta's dominance of the 
spectrum occurred at less than 0.5 keV where interstellar absorption makes 
observations difficult. The amount of emission from the ejecta was comparable 
to the shocked ISM at ^ 1 keV, but as with discussion on the iron abundance, 
varying the value of n, which varied the contribution to the a. 1 keV emission 
from the shocked ISM, could compensate for any changes in the ejecta emission. 

In low resolution we were not able to decouple the components. 

Using much the same argument for providing a lower limit to the distance we 
can put an upper limit on the size of 6 of 'v 40. Any value higher than that would 
produce shock velocities that were too low to produce a 5 keV plasma during the 
lifetime of the remnant. 

Blast Energy 

Knowledge of the blast energy has a direct application to accessing whether 
the progenitor of MSH 14-63 is a Type I or Type II supernova. The issue of the 
association of MSK 14-63 with the OB star group also relate, the question of the 
type of SN explosion that formed the remnant. The association would strongly 
suggest that the MSH 14-63 was a Type II supernovae (Shklovsky, 1962), since 
Type II supernova are thought to occur in massive stars. 

The results of the standard NIE model with 8 fixed at unity which are given 
in Table 2 would indicate that the blast energy, E q , was £ 6 x 10^ ergs. How- 
ever, as the discussion in Chapter 2 pointed out, E q is a very sensitive function 
of n- Fixing the abundances to solar, for example, moved the inferred value of 

51 

the blast energy up to 'v 1.2 x 10 ergs. This large range of uncertainty in E q 
along with the large range of uncertainty in the mass of the ejecta, M g does not 
allow an unambiguous assertion as to the type of supernova explosion. 



[Fe XIV] Emission 

Lucke et al. (1979) have made measurements of the [Fe XIV]X5303 emission 

-14 -2 -1 

line in RCW86. They observed a flux of 2.28 + .98 x 10 ergs cm s from 

summing four 5 arcmin circular fields in the SW part of the remnant at the shock 

-12 -2 -1 

front. The dereddened flux was given as 1.4 x 10 ergs cm s . In LIE, [Fe 
XIV] emission peaks rather sharply at 2.1 x 10^ K. For CIE plasmas with 
T g ^ L5 x 10^ K or £ 3.0 x 10^ K the emissivity of the X5303 line drops 
precipitously as a function of temperature. The authors- interpreted their ob- 
servation in terms of the two-temperature model fit to the x-ray data. They 
ignored the high-temperature blast-wave component, because under the assumption 
of CIE the amount of [Fe XIV] emission would be negligible at the measured tempera- 
ture of 6 x 10 7 K. Rather, they suggested associating their flux with the alleged 
low-temperature x-ray component (2.5 < T^ < 6) or an even lower temperature 
component which they argue would be difficult for x-ray detectors to detect 
because the bulk of this emission would be attenuated by photoelectric absorption 
in the intervening ISM. 

Though it is possible that these low temperature components exist, we would 
like to put forth an argument that ignoring the high-temperature blast wave is 
not justified. Because of the time-dependent nature of the ionization behind the 
shock, a shell of [Fe XIV] is produced which has sufficient emission volume as to 
produce a flux comparable to that observed. 

The unattenuated [Fe XIV] X5303 flux observed at the earth, S, from the blast 
wave of a remnant at a given time is 


4n/ S n 2 P,.(T )A_ X, . x 2 dr 
J e 14 e re 14 
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where R is the shock radius, n is the electron density, P., is the emissivity of 

the line assuming a solar abundance for the iron (from Allen, 1976) and 100% 

concentration of the iron in the (Fe XIVJ ionic state, Aj, g is the abundance of 

iron in units of solar, is the fractional concentration of the [Fe XIV] ionic 

state, r is an integration variable, D is the distance to the remnant, and f^ is 

2 

a correction for a beam size less than the size of the remnant. X,. and n are 

14 e 

the most rapidly changing quantities behind the shock so we approximate the 
equation for S as 


S 


: ,v f „b W'-’V 


where <X^> is the emission weighted average of the X^ and is given as 

X A A 


<X 14 > 


C n e x !4 r dr 

/ 8 n 2 r 2 dr 
o e 


The [Fe XIV] X5303 line results from the forbidden transition 
2 2 2 2 

3s 3P ~ P 3/2* Population of the P^ 2 level results from direct electron 

excitation of the ground state or from cascades following electron excitation to 
higher energy levels. In general, the emissivity, for an atomic line from an 
element with nuclear charge Z and ionic charge z, in an optically thin plasma of 
temperature T , and resulting from electron collisional excitation is given by 



where X„ is fractional ionic concentration, (n_/n u ) is the number fraction of 

l* f 2 tx 

the element relative to hydrogen, (E/I^) is the energy of the emitted photon in 

rydbergs, - is the effective collision strength, E is the excitation energy, 
ert ex 
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and Tg is the electron temperature in units of 10° K. To arrive at P^, the frac- 
tional ionic concentration is set to unity, the exponential factor is set to unity 

because kT >> E , and the elemental abundance taken from Allen (197C). A 
e ex’ 

value for of ^3 is obtained by extrapolating from a graph given in Nussbau- 

mer and Osterbrock (1970). ^ e ff Includes cascades following excitation to 
higher levels. The resulting expression for is 


P 17 « 3.51 x 10 
14 


-24 T e ” 1/2 Z eff -1 3 

— r — ergs s cm 


io 6 k 


We will also include X5303 photons produced following ionizations of Fe XIII 
to FE XIV. Using branching ratios given by Nussbaumev and Osterbrook (1970), we 
calculate that 46% of the ionizations will result in X5303. Taking ionization 
rates from Raymond and Smith (1977), we can calculate an analogous expression 
to P u for P 13 : 


P 13^ T e ^ * 3.4 x IQ' 26 ergs s 1 cm 3 . 

7 13 -5 

We use a nominal value of 5 x 10 for T and 1.4 x 10 cm for C . For the 

e ev 

calculation of <X^> we will assume that at the shock, at least, a low n model as 

valid. For n ■ 0.16 and x » 1.2, <X^> 0.18 and <Xj^> 0.16. Wt estimate 

f , , which is fraction of emission volume observed, to be 0.07. This calculation 
ob 

of f ^ incorporated effects of limb brightening. These values given S equal to 

n m-14 . -3 -1 

9 x 10 ergs cm s 

This number is less than the measured value by a factor of ^13, but three 
good possibilities that could bring the calculation more into agreement are: 

1) The iron abundance of the shocked material could be high by a factor of 2-3 


times solar. 



2) Lucke et al. observed in the region of the most intense x-ray emission* This 

implies that the t jst-shock densities might be higher in this region* Higher 

densities shrink the X^ by a factor of n \ but they increase the surface 

2 

brightness, which is proportional to C^f ^ by n^ • The resulting value of S 

scales like n • 
e 

3) The Lucke et al. measured value for S is dependent on dereddening. Th.y 
took 2 mag of extinction from che measurements by Westerlund (1967) which 
assume thac the remnant is In the OB Association. If it is closer as the 
x-ray data suggest, the reddening correction would be lower whan the one used. 

Given the uncertainties in the input parameters, statistical uncertainties, 
and uncertainties in the dereddening, this predicted flux might be consistent with 
the one measured. At the very least, the [Fe XIV] flux from the high temperature 
blast wave component cannot be dismissed as lighuly a g Lucke et al. suggest if 
the plasma is not in CIE. For comparison, we note that the value of S assuming 
CIE in a 'v 5 keV plasma is decreased ^ 20 orders of magnitude from the NIE result. 

Conclusion 

In summary, the results of this chapter are the following: 

1. A significant non-Coulomb, ion-electron interaction is occurring in the rem- 
nant, presumably at the shock front. 

2. The data can be fit by a model with little or no emission from any source 
other than the shocked ISM with an age which Is consistent with that of MSH 
14-63. 

3. Anomalous abundances of heavy elements and possible discrepancies in the 
centroid and shape of the Fe feature could be explained by inhomogenities 
in the ISM density or by ejecta that have come to thermal equilibrium with 
the shocked ISM. 

4. Under the assumptions of the model, the distance to the MSH 14-63 is incon- 


sistent with distance measured to ar OB association in the direction of the 



remnant; however, plausible deviations In the model assumptions might explain 
the discrepancy. 

Future observations with good combined spectral and spatial resolution, 
especially at Fe would be helpful in addressing the issues raised here. 
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Table 2 

N1E Model (T. = Tj) 


best fit 

~ limits 

C„ (10 19 cm* 9 ) 

1.4 

1.6 

1.2 

ti (foe l/9 cm* 8 ) 

0.16 

0.12 

022 

fi (Af 0 foe* 1 ) 

1 (fixed) 

- 

- 

T (kiloyears cm* 1 foe _l/9 ) 

1.2 

1.3 

1.1 

(10 81 cm* 8 ) 

4.4 

4.7 

4.1 

(solar) 

0.7 

! 

0.6 

0.B 

As (solar) 

3.4 

4.1 

2.B 

Xp e (solar) 

39 

4.6 

2.6 

X* (v = 53) 

1.3 

1.4 

1.4 

a (foe 179 cm) 

0.76 

0.33 

1.7 

Eq (foe) 

0.09 

0.01 

0.6 

n 0 (cm* 9 ) 

0.2 

0.4 

ESI 

fa (10 s years) 

0.9 

0.4 

1.9 

D (kpc) 

0.7 

0.4 

1.6 


cb. foe a 10 s1 ergs (fifty-one ergs; 9 Heins Bstbe) 
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Table 3 


NIE Model with K° lines 


best fit 

r w (10‘ 3 cm' 6 ) 

1 2 

17 (foe l/8 cm"*) 

0 30 

^(Af e foe-») 

1 (Axed) 

r (kiloyears cm'* foe~* /s ) 

1.1 

JVj, (10** cm"*) 

3.5 

X* (* = 53) 

1.3 

a (foe 1/5 cm) 

2.6 

£9 (f°«) 

2.4 

n 0 (cm* 8 ) 

0 11 

t s (10 9 years) 

3.1 

D (kpo) 

2.9 


Figure Captions 


Figure 1 

(a) shows a fit to the data using a two-termperature CIE model. Histogram 
represents best fit model and the crosses are the data points, (b) shows only 
the high temperature CIE model. 

Figure 2 

An illustrative example of electron temperature profile when only Coulomb heat- 
ing is considered is given in (c) from calculations by NIE model. Cox and Ander- 
son (1980) approximate analytic solution is denoted by the broken line. The 
density and mean temperature profiles are shown for comparison in (a) and (b) , 
respectively. Histograms are calculations from the NIE model and the smooth 
curves is calculated from the analytic solution for a adiabatic remnant. 

Figure 3 

The best fit model (histogram) using NIE model is shown with the data points. 
Insert is the projection of X^ space on the n - t plane. Contours show the 
change of from the best fit. Both 4.6 and 10 are shown. 

Figure 4 

Same as Figure 3 but now lines have been inserted at the Ka energies of S, Ar, 
and Fe (2.4, 3.0 and, 6.7 keV, respectively). The abundances of the elements in 
the NIE model were set to solar. 


Figure 5 

These graphs represent fits to the data at photon energies £4.5 keV where the 
non-iron contributions to the model have been subtracted to emphasize the fit 
of the model at the Fe Koi feature, (a) shows the best fit when n is fixed at 0.2, 
(b) is the best fit letting all parameters vary, and (c) is a fit with a CIE model. 
Fitted abundances, Ap^ are given in the figure. 
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4*3 The Observation of RCW 103 


The first report of x-ray emission from RCW 103 was made by Tuohy st al. 

(1979) using data from the HEAO A-2 low energy detectors (LED). The detection was 
made at photon energies ^ 3 keV. Fits to the pulse height data using a collisional 
ionization equilibrium (CIE) plasma model allowed temperatures in the range of 
1.3 - 2.0 x 10^ K. Results from tiu. HEAO A-2 medium energy detectors (MED) 
were ambiguous because of source confusion. An upper limit to the 2-10 keV 
flux was 1 x 10 ^ ergs cm ^ sec As Tuohy et al. point out, this upper limit 
does not exclude any temperature that is consistent with the low energy spectral 
data. 

X-ray images have been obtained from the Einstein Observatory. Data from 
the High Resolution Imaging (HRI) experiment (Tuohy and Garmire, 1980) and the 
Imaging Proportional Counter (IPC) experiment (Tuohy, private communication) 
exhibit a limb brightened shell. The surface brightness around the ring is not 
uniform - the most intense region is in the southeastern part of the remnant. 

The angular radius of the x-ray emission is ^9 arcminutes. Radio maps of the 
region also show a ring-like structure of comparable radius (Goss anc Shaver, 

1970). Optically, filaments exist in the southeast northwest parts of the 
remnant. They are contained within a projected radius which is smaller than 
the radius of the x-ray/radio emission. 

An interesting feature of RCW 1C3 Ls the presence of a x-ray point source in 
the geometric center of the remnant (Tuohy an'. Garmire, 1980). No detectable 
radio emission has been observed from this point source (Tuohy et al. , 1983) 
suggesting that it may be of a different character from the pulsars that exist at 
the center of some other supernova ^emnants (SNR). Tuohy and Garmire suggest that 
it may be the first detection of thermal radiation from a cooling neutron star. 

A distance to RCW 103 of 3.3 kpc is estimated based on the kinematic models 
of the neutral hydrogen in the spiral arms and absorption profile of RCW 103* s 



radio continuum at 21 cm (Caswell et al. , 1976). This distance estimate along 

with the measured angular diameter of the x-ray emission would determine the shock 

radius to be 4.5 pc. This size would indicate that RCW 103 is a young or possibly 

middle aged SNR. Assuming the kinematics of the shock are described by a Sedov 

1/2 

blast wave, the remnant size would imply an age of 'v 600(n^/E^) years, where 

-3 

the vde density of the interstellar medium (ISM), n^, is in units of 1 cm , 
and the supernova blast energy, E^, is in units of 10^* ergs. 

Westerlund (1969b) has suggested that the progenitor of RCW 103 might have 
been a member of an OB association 3.9 kpc away that lies along the line of sight 
to the remnant. The association with an OB star group \ 'uld also suggest that 
the supernova was a Type II event since Type II supernovae are thought to be due 
to massive, early type stars (e.g. Shklovsky, 1962). 

In the previous chapter we presented an analysis of the HEAO A-2 data of 
SNR, MSH 14-63 (RCW 86), using a model developed to account for non-equilibrium 
effects in SNR (NIE model). In this chapter we use the model to fit the 
RCW 103 data collected by the Solid State Spectrometer (SSS) experiment on the 
Einstein satellite. Because the SSS has good spectral resolution below % 4 keV 
(AE/E -v 10%), line features can be more readily distinguished with the SSS than 
with proportional counter detectors flown previously. These line features, each 
associated with a blend of lines from a single element, are seen prominently in 
the spectra of most young SNR's. In the SSS, K a lines from Mg, Si, and S are the 
three most striking features. Lines from n*3 to n*2 transitions in Fe also 
produces a large amount of the flux at around 1 keV. io properly infer elemental 
abundances from the intensities of these features, a proper accounting of the 
ionic concentrations must be made. Previous work modeling the RCW 103 SSS data 
have used a CIE plasma model (Becker, private communication). Two plasmas et 
different temperatures were required to explain the data. In addition, the CIE 



model required chat the abundances , relative to solar, for the elements of Mg, 

Si, S, and Fe be \ 2. Because the plasma may not be in CIE the calculated ionic 
concentrations and hence the inferred elemental abuidances may be in error . 

Models which assume non- ionization equilibrium (NIE) in a solar abundance, shod* 
heated, plasma produce a spectrum, that when convolved with a detector response 
function, can emulate a second, lower temperature component and produce anomalous 
line strengths relative to a CIE model (ltoh, 1977, 1978, and 1979; Nugent, 1982; 
Shull, 1982). 

The reflection efficiency of the telescope on the HEAO A-2 observatory drops 
sharply for photons above ^ 4 keV. Consequently no high energy data were obtained 
with the aforementioned imaging observations (IPC/HRI) or during the SSS 
observations discussed above. This ignorance of the high energy spectrum has 
serious consequences for the analysis performed here. Because line emission from 
heavy elements is the dominant emission mechanism in the part of the spectrum 
where the SSS is sensitive, interpreting the spectral data requires a complex 
model; however, above v 4 keV the continuum dominates, so measurement of the 
spectrum at these energies is a less model-dependent measure of the electron 
temperature and reduced emission volume of the shock heated ISM. Addressing the 
issue of electron-ion equilibration behind the shock becomes more difficult as 
result of this lack of dati. Measurement of the continuum in this energy range 
has been definitive in demonstrating that some non-Coulomb electron-ion interactions 
are present in Cas A (Pravdo and Smith, 1979), Tycho (Pravdo and Smith, 1979), 
and MSH 14-63 (Nugent et al., 1982). 

Despite a growing body of evidence that these non-Coulomb processes exist 
at the collisionless shock fronts of young supernova remnants, there still does 
not exist a detailed model of the process, so that we are ignorant of how the 
magnitude of the effect might vary over a plausible range of SNR shock parameters. 



It is known, for example, in the collisionless shocks which occur when the solar 
wind impinges on the magnetosphere of the earth, that the electron temperature 
is lower than the ion temperature* Even though we have a bias towards models 
that have electron-ion equilibration, we will search for evidence of this in the 
SSS data on RCW 103* 


Observations 

The spectral data were obtained from observations with the Solid State 

Spectrometer (SSS) experiment on the HEAO-2 spacecraft. The SSS experiment is 

described in detail in Joyce et al. (1978). Briefly, the experiment has an 

2 

effective collecting area of 200 cm and a circular field of view with a diameter 
of 6 arcmin. Spectral data are collected in the photon energy range from ^0.5 to 
^ 4 keV. 

Because the field of view of the spectrometer is smaller than the angular 
size the remnant (9 arcmin), a number of different observations are required to 
obtain the spectrum from the entire remnant. Specifically four observations 
were performed. Figure 1 is an x-ray surface brightness map showing the positions 
and approximate extent of the field of view of the different observations* For 
each of the observations. Table 1 gives the co-ordinates of the center of the 
field of view, the epoch, and the integration time. The uncertainties in the 
positions are £ 1 arcmin (Symkowiak, private communication). 


Comparison of Observations 

The fields of view of the four observations of the RCW 103 remnant do overlap 
to some extent, but they all sample a different part of the remnant. This gives 
some limited spatial resolution to the spectra. Before fitting the data we shall 
statistically test whether the four data sets are consistent with a single, general- 
ized spectral model. By a generalized model, we mean a set of numbers, M ^ , where 
the subscript refers to the jth pulse height channel. Because the surface bright- 



ness varies over the remnant and because the amount of the remnant subtended by the 

field of view during a certain observation varies we will allow the observations 

to have overall scaling factor, s^. The subscript in refers to the ith 

observation. 

2 

The X statistic used can be written as 


X - 


I 

V 

4* 

i-i 


j 

z 

J-i 


(s i M r D i/ 

2 

°1J 


where is the count rate in the jth pulse height channel of the ith observation* 
is the lo uncertainty in D^* J Is the number of pulse height channels and 
I is the number of observations. 

The analysis is complicated by a number of effects. First* because certain 
emitting regions were common to two o.: more observations* the data sets were not 
completely independent. For simplicity, though* it was assumed that they were. 

A second problem concerned the thickness of an absorbing layer of ice that 
accumulates on the window of the SSS detector. The thickness* hence the amount 
of absorption of the impinging x-rays varied from one observation to another. The 
correction of the data for this absorption* in principle* is dependent on the 
c’noice of model for the input spectrum. Since no spectrum was assumed for this 
statistical test, an approximate correction was made by dividing the data for 
each pulse channel by the absorption factor for photons with energy corresponding 
to the nominal energy for that channel. A final problem with the data concerns a 
background component that is known to vary. This problem is most severe in 
channels which correspond to photon energies v 2.5 keV. We felt that the shape of 
the background contamination in these channels is different from the nominal one 
that has been calculated (private communication* Svmkowiak) . More attention will 
be given to this issue later in this paper, but for this comparison her ? will 
ignore the data above 2.5 keV. 



2 

The result of the analysis is displayed in Figure 2. The reduced X for the 

2 s 

fit was 1,28 for 117 degrees of freedom* The chance of being ^ 1*28 for 117 

degrees of freedom if .he model were correct ia 2*3X. Under the usually criteria 

2 

of "90i confidence" the model falls the X test and does suggest the possibility 

of intrinsic differences of the spectrum arising from different regions in the 

2 

RCW 103 remnant. The slightly high X could be also due to the systematic errors 
outlined above. Therefore, for simplicity, we will adopt the latter assumption for 
the remainder of this paper. Given that the spectra are very similar, it isques- 
tionable whether the exercise of modeling each data set separately would produce 
any meaningful differences. 

As was mentioned in the introduction* x-ray images of RCW 103 reveal that there 

exist surface brightness differences over the surface of the remnant* Because of 

uncertainties in the precise position of the SSS fields of view and the lack of 

any detailed data on the x-ray morphology aside from a plot, we cannot confirm 

that the flux received in the SSS bandpass is consistent with the measured surface 

brightness with the IPC or HRI detectors. All three Instruments have different 

bandpasses - the SSS being the smallest, so if the SSS and IPC or HRI measurements 

of the relative surface brightness were inconsistent that would imply that spectral 

differences outside of the SSS bandpass were responsible for the differences in 

remnant’s x-ray morphology. Alternatively, if that SSS flux is consistent with 

2 

l PC /HR I result, then variations in the emission integral, [/n *“dll. would be 
indicated. Lack of variation in the SSS spectrum implies that the ionization time, 
,’n^dt, is constant at different positions around the remnant which would further 
implv that there are no density variations and that variations in filling factors 
may be responsible for surface brightness variations. 

In order to aid future workers on the x-ray morphology of RCW 103, we give 
the relative do iced flux of the four observations in Table 2. The fluxes are 


normalized to the flux of observation f*l. 
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Review of Parameters 

Before discussing results of the model, we briefly review for the reader the 
definitions of the parameters used. The input parameters to the model and the 
units we will always use are the initial blast energy, E 0 (units of 10^ ergs), 
the density of the interstellar medium around the progenitor star, n Q (units of 1 

_3 

cm ), the amount of ejected mass from the progenitor star, M g (units of Mg), the 
age of the remnant, (units of 1000 years), the abundance, relative to hydrogen, 
of elements in the ejecta (normalized to the measured solar value), the abundance, 
relative to hydrogen, of elements in the ISM (normalized to the measured solar 
value) , and the distance to the remnant , D (units of kpc) . These parameters are 
mapped to a more convenient set of parameters for purposes of the doing spectral 
fitting. There is a scaling parameter defined as 


a 


( 



1/3 


and a number of "reduced" parameters: 


n * (E n 2)*^ = a = an , 
o o o o’ 


M 


-1-2.. 


8 = f- * n ~a X * 


and 


T = 




oTV 


A fourth reduced parameter is the reduced emission volume, C , which is defined 
as 
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c 

ev 


r n 2 dV 


n, T, B, and C ev are derived from fits to the spectral data, a is then 
determined from the average x-ray surface brightness by 

„ _ £'(n,B,T)» 2 f 
C 

ev 

where is determined from n» 6. and t and f is the fraction of solid angle of the 
remnant filled by shocked plasma. The value of f will nominally be asumed to have 
a value of unity. With a, the input parameters are given by 


2 2 -1 
E - on M »ctn6 t_*aT n*un 
o e Jo 


and the distance is given by 


D 


R g * (B» t)o 

§ 


where R ' is a constant determined from the fitted value of B and t. 
s 

When t £ 16B 5 ^ 6 the remnant is in a Sedov phase and the shock temperature, 
T g , can be given as a function of x by 


3.2 t- 6/5 keV 


Fits to Spectral Models 

As was discussed in the section concerning comparison of observations, there 
are systematic problems with the background data above 3.0 keV. In the case of 
the RCW 103 data there were problems reconciling the data with the fitting models 



at photon energies £2.5 keV. We will first present results, considering fits 

to the data below 2.5 keV and then discuss the residuals above 2.5 keV. 

We will first discuss models similar to those used by Hamilton et al. (1982). 

As discussed above these models use Sedov hydrogynamics and are equivalent to 

this model only in the "far Sedov" limit (t >> x g ). The models were fit with 8 

fixed at 2 (8 must be finite for numerical reasons.) rad no ejecta emission was 

considered. Models where the ejecta were included are discussed below. The 

results of the fit with the "hot" electron (T = T ) model are given in Table 3 

e m 

and the dtat with the model are shown in Figure 3b. Likewise, the "cold 1 * electron 

(T ^ T.) model results are given in Table 3 and the data with the model are 
e i 

displayed in Figure 3a. For both of these fits t/t £ 12, thus assuring our desire 

s 

to be in a "far Sedov" limit. The parameters that were varied were C fiv , n* t, 
and the abundances for magnesium (A^) , silicon (Ag^), sulphur (A g ), and iron 
(Ap e >. The rest of the elemental abundances were fixed at their solar value. 

On the assumption that the emission is from a uniform ISM shocked by a 
Sedov blast wave, we have inferred the input parameters by the method described 
in the last section. In addition to the filling factor, f, there is a related 
effect that the SSS field of view is smaller (6 arcminute) than the extent of 
RCW 103 (9 arcminute), and, thus the measured C ^ given in the tables is an under- 
estimate. Most of the exposure of the data comes from observations 3 and 4 where 

the remnant filled the entire field of view, so, as an approximation, we corrected 
C^ v before calculating a by the ratio of the solid angle of the remnant to the 
solid angle of the SSS beam. We took that ratio to be ^2. The value for a and 

the inferred input parameters are given in Table 3. 

2 

The fits to the data are not of high quality. The reduced X obtained would 

2 

result < 0.1% by chance. It is possible that the high X are in part due to 
improper modeling of detector output to a given input spectrum indlucing back- 



ground effects, but, even if all the problems were due to detector modeling as 

opposed to source modeling, there is no reason to believe that our best fit 

parameters are close (compared to the statistical uncertainties) to the best fit 

values given ideal detector modeling. Despite this caveat, we performed on the 

2 

"hot" electron model a contraction of X space onto the n-x plane. Every point 

2 

on this plane represents the best fit X if n and x are fixed at the values given 

at that point and all others are allowed to vary. The results are shown in the 

2 

insert of Figure 3b. Two contours are given. The first is for values of X which 

2 

exceed the minimum X by 4.6. This corresponds to 90% confidence region for 2 

2 

interesting parameters (Lampton et al., 1976). Another contour at AX * 10 is 
given. Two representative points near the extremes of his region are also shown 
in Table 3. These illustrate rough statistical limits for the various parameters. 

We note that the "valley" of X 2 on the n-x plane roughly runs along a line of 
constant effective ionization t i * , nx. The fits are snesitive to the shape of the 
Ka lines of the principal elements. Each K a blend is composed of emission from 
the hydrogen-like ion, helium-like ion, lithium-like ion, etc. Each ion emits at 
slightly different energies that are roughly ordered by the ionic charge. The SSS 
does not, of course have the capability to resolve these blands of lines, but they 
can be inferred from the shape and centroid of the line which, assuming kT g is 
D'urh greater than the energy spacing, will be dependent on the relative concen- 
tration of the ionic species. These ionic concentrations depend somewhat on the 
temperature history, but for the temperatures and elements of interest here, the 

ionic structure is a more sensitive function of / n dt. This is the reason that 

e 

2 

X contours remain close to iso-nx contours. 

We will first discuss the significance of the different parameters individually. 
Following that, we will make some general comments and summarize the main points. 



Contribution from the ejecta 

The possibility of an ejecta emission component resulting from a "reverse" 
shock wave has been discussed by McKee (1974) and Gull (1973, 1975), As has 
been discussed this reverse shock component manifests itself in one-dimensional 
hydrodynamic codes like the one used in this model, so, in principle, we can use 
our NIE to search for the evidence of this component in the x-ray spectral data* 
When reverse shock emission from the ejecta was included in the model, best 

2 

X were obtained when the ejecta did not contribute significantly to the 
measured flux. Either the ejecta had a small fraction f the emission volume or 
the ejecta component had a temperature £ 100 ev-a temperature to which the SSS 
is insensitive. It is possible, however, to find models that do have a signifi- 
cant ejecta component that are not in gross disagreement with the data. To 
demonstrate this a search was made in a restricted part of parameter space for 
which t < 4T g . 6 was allowed to vary, along with the parameters mentioned above. 
Our spatial and spectral resolution were not sufficient to deconvolve the ejecta 
contribution to line strengths from the ISM contribution. Hence the abundance 
parameter for a given element in the ejecta and the abundance parameter for the 
same element in the ISM are highly correlated in the fitting process. One can 

be arbitrarily decreased and, to compensate, the other can be increased without 

2 

significantly changing X . This problem can be avoided by fixing one of the 

abundance sets. In most modeling runs the contribution of the ejecta turned out 

to be a smaller fraction of the emission, so the ejecta was chosen to be the 

fixed set. Solar abundances for the ejecta were used; however, this might be 

a poor choice, especially if the progenitor was a Type II supernova. 

The results of this search for "hot" and "cold" electron models are shown in 

Table 5 and Figure 4. In the case of the "cold" electron model a local minimum 
2 

X 4 * in parameter space was found, but for a "hot" electron model x needed to be 



constrained in order to satisfy the condition of being < 4t . Both fits give some 

3 

2 

what worse but comparable X to models where the ejecta was ignored. 


Electron heating at shock fronts 

The issue of whether electrons are heated by non-Coulomb processes in the 
chock of RCW 103 is not clearly resolved by the fits to the SSS data. One reason, 
mentioned before, for this ambiguity is the inability to directly measure the 
high energy continuum which gives a relatively model-independent measure of the 
electron temperature. That inability is related to the insensitivity of the SSS 
to photons with energies ^ 4 keV and is a general problem when relying solely on 
data from the SSS experiment. Assuming the results of the NIE model are roughly 
correct there is a more fundamental problem, intrinsic to the RCW 103 SNR. 

RCW 103 is apparently in a part of n-B-i space where the "cold" and "hot" models 
are beginning to become indistinguishable. The age, t , at which the electron 
temperatures predicted by a "cold" electron model and a "hot" electron model 
begin to become indistinguishable is given by Cox and Anderson (1982) as 


t 

eq 


3000 E 5 


3/14 



years 


where is the balst energy in units of 10^* 

-3 

sity in units of cm . The ratio of the age, 
in terms of our parameters, is given by 


ergs and n is the ISM particle den 
o 

t, to t is independent of a and, 


t t 5/14 
— - 3 n 
eq 

For n • 0.5 and t ■ 3.6, t/t = 0.94. The fact that the best fit points for both 

eq 

models lie roughly in the same region of parameter space is consistent with the 



fact that the models a . beginning to converge in that region. For t >> t the 
differences in the two model spectra become negligible, but at t ^ t only the 
differences in the predicted electron temperature and hance the continuum spectra 
becomes negligible. The fact that the fits are not completely similar can be 
explained by the different electron temperature histories of the shock gas. In 
the "hot" electron model, the election temperature of a particular gas element 
has been decreasing monotonically ever since that gas element has been shocked, 
while during some or all of that time the "cold" electron model predicts that the 
temperature has been increasing. Since the evolution of the ionisation structure 
is dependent on the evolution of the electron temperature, it is reasonable to 
expect differences between the two best fit points in parameter space, 
particularly the elemental abundances. 

Elemental <*Dundances 

The fitted ISM elemental abundances in the "hot" electron model are consistent 
with a solar set. This differs from the results using CIE models by Becker 
(private communication). The abundances required by the CIE models are a factor 
of 2-3 (Table 5) higher relative to solar. The abundances determined by the 
"cold" electron model deviate from solar by up to factor of % 4. 

Since a solar abundance set is nominally expected from a interstellar gas, the 
results of the NIE "hot" electron model are pleasing. The abundance results may 
be closer to solar than are shown in Table 3. For historical reasons an Allen 
(1976) abundance set was used. More recent work by Meyer (1979) shows some 
significant deviations from Allen for some elements. The conversion between to 
the two abundance sets is complicated by the fact that the continuum flux is 
dependent on the elemental abundances, particularly C, N, and 0; however, as an 
approximation we will correct each element by its ratio in the two abundance 



sets. With that correction, the best fit abundances for Mg, Si, and Fe are all 
within 5% of solar. This makes the Meyer solar value consistent for each of these 
elements. The only statistically significant deviation is for S which has Meyer 
abundance of 1.2. 

Not all results from applying N1E models to SSS data of young SNR's have 
demonstrated solar abundances. Work on Tycho’s SNR (Shull, 1982), for example, 

has shown overabundances still are predicted by NIE models. The distinction may 

* 

result from the dominance of the role piayed by the ejecta in the x-ray emission, 
since the ejecta is thought to be metal rich. In addition to the emission from 
the ejecta following heating from a reverse shock, ejecta may also ix with hot 
ISM and heat by conduction. 

Distance, density, and interstellar absorption 
The distance measurement is consistent with the 3.3 kpc estimate made by 
Caswell (1976). It also is consistent with the distance measured to Westerlund’s 
OB star group, strengthening the argument for association of the progenitor of 
RCW 103 with that group. 

The height of RCW 103 off of the galactic plane is 23 pc assuming a distance 

-3 

of 3.3 kpc. The inferred particle density of ^ 1 cm is reasonable for that 
position. 

Finally, assuming a distance of 3.3 kpc, the measure of N^, would imply a 

-3 

mean density of ^ 0.4 cm in the Intervening ISM. Again, this is a reasonable 
value. 


Age 

The age ot RCW 103 has consequences for the issue of the lack of historical 
record of the supernova which formed RCW 103 and of the cooling of central point 
source seen by Tuohy and Tarmire (1979). We will adopt for this discussion the 



age, with uncertainties, determined from the "hot" electron model. First, to put 
this work in perspective, we note that limits can be put on the age of RCW 103 
without resorting to the use of the elaborate NIE model. Taking the radius of 
the remnant to be 4-5 pc and assuming that the average velocity f the ejected 
material is £ 0.1 c would put a lower limit to the age of * 140 years. From 
crude spectral modeling of x-ray data a lower limit for the shock temperature 
can be determined to be ^ 0.2 keV. That limit on the temperature along with the 
measurement of the size of the remnant can be used to assert an upper limit to 
the age of ^ 6000 years. 

Due to its apparently young age and close proximity to the sun, RCW 103 
could be expected to be a historical supernova. As pointed out by Clark and 
Stephenson (1976) if RCW 103 is younger than 600 years the initial supernova 
event could have escaped attention by Chinese astronomers, since the capital of 
China was at Peking for most of the time after AD1280. Peking is too far north 
(latitude ^ +40) to able to see RCW 103. Chinese records become more scarce 
before ^ 200 BC, so, if the age of RCW 103 is £ 2200 years, the non-detection of 
its supernova, could again be explained. At intermediate times, records existed, 
the capital was farther south, and the effect of the earth’s precession made 
objects in the direction of RCW 103 more visible in northern latitudes. From the 
NIE model, RCW 103 is most likely of this intermediate age. To emphasize the 
point that the Chinese were able to monitor objects in the direction of RCW 103, 
they did record supemovae in AD 185 at a { 15*\ 6 % -62° (1950) and in AD 393 
at u £ 17* 1 , 6 % -40° (1950). They should have had an unobstructed view of the 
RCW 103 supernova, which is at a ^ 16 , 6 % -51° (1950), if the event took place 
at the tim'* determined by the NIE model. 

As also pointed out by Clark and Stephenson (1976) even if the RCW 103 
supernova took place during a time when Chinese could observe it, the likelihood 



of detection could still have been reduced if RCW 103 was near conjunction with 
the sun at the time of the event. We also point out that Cas A, which is thought 
to be at the same distance as RCW 103, most likely escaped detection as a his- 
torical supernova. If anything, Cas A was a more likely candidate for detection 
since ii was probably more recent and definitely at a much more favorable 
declination for viewing in the northern hemisphere. 

A second issue related to the age of RCW 103 is the cooling rate of the point 
source in the center or the remnant. Under the assumption that the object is a 
neutron star and that no mechanism exists for renewing the internal energy of 
the core, then knowing the age of the neutron star could be a test for different 
suggested cooling mechanisms (Bahcall and Wolf, 196S; Glen and Sutherland, 1980; 
Normoto and Tsuruta, 1981; Richardson et al., 1982). The NIE model can only be 
an effective discriminate of cooling mechanisms if the cooling time scale is 
thin this allowed range for the age of RCW 103. This probably is not the 
case. Neutrino cooling following pion beta-decay is the fastest cooling mechanism 
with cooling time scales ^ 1 year. This can be easily ruled out for RCW 103's 
neutron star. Recent calculations by Richardson et al. which do not assume a 
pion concentrate in the core give cooling curves that are consistent with the 
age measured here and the effective temperature of 2.2 x 10^ K measured by Tuohy 
and Garmire. In principle, the age and effective temperature could be used 
within the framework of models - such as Richardson et al. - to determine the mass 
of the neuron star; however, the uncertainties in the measured quantities were 
too lar^e to allow this. 


Blast Energy 

The value determined for the blast energy. 


v* 6 x 10 


50 


ergs. 


is somewhat low 


compared to the numbers that are usually associated with supernova blasts 



(^0.4 -4 x 10^* ergs). Some care must be taken before accepting this upper limit, 
however. Because the inferred value of the blast energy scales like it would 

only take an increase of *v 10% in the fitted upper limit of n to account for the 
discrepancy. 

Since Type I supernovae are thought to have lower blast energies (E^j *v 0.5) 
then to Type T I (E^ % 1) , the low fitted value of would be weak evidence for 
the progenitor of RCW 103 being a Type I supernova. 

Summary 

Our inability to settle the question of electron-ion equilibration at the shock 
front may be due, as suggested, to intrinsic properties of the RCW 103 remnant; 
however, this ambiguity and the ambiguity on the question of reverse shocked ejecta 
are exacerbated by poor knowledge of the spectrum above 2.5 keV. This part of the 
spectrum is dominated by continuum emission and as a result gives a measure of the 
electron temperature independent of many of the systematic problems of the NIE 
model. It also gives a direct measure of the reduced emission measure of this high 
temperature component which provides a baseline for interpreting the "excess" of 
flux at 1-2 keV. As can be seen from the analysis of the RCW 86 data, this infor- 
mation constrains the permitted volume of parameter space severely even without 
detailed knowledge of the line emission. 

Knowledge of this component could, for instance, help to distinguish between 

"cold" and "hot" electron models. The "hot" model would predict a shock temperature 

13 -5 

or 'v 1-2 keV and reduced emission volume of ^ 1 x 10 cm while the "cold" model 
woulc have the high temperature be 'v 0.7 keV and the reduced emission volume be 
^ 3 x 10^** cm 

Am important future observation will be the determination of the behavior of 
the spectrum at photon energies above 2.5 keV. Source confusion has thwarted 
efforts to do so with HEA0-1 and the high energy cutoff of the x-ray mirrors on 
HEA0-2 prevented the SSS from operating in that region. AXAF, BBXRT, as well as 



more modest shuttle and Explorer missions should have the capability to make this 
measurement . 

Despite the poor coverage of the spectrum, the SSS data on RCW 103 does, of 
course, provide a good opportunity to study the intensity of the line emission 
and, hence, deduce elemental abundances. On this note we consider the NIE model 
to have been a success. Assuming the "hot" electron model is correct, the data 
dc support the idea that the x-ray emission is arising from shocked ISM with 
elemental abundances roughly the same as solar. The f hot' electron model is 
preferred because the "cold" electron model has been ruled out in other yound SNR's. 

Finally, we note that there are indications that RCW 103 may have up to four 
unique or almost unique characteristics relative to other young SNR's. Wich 
parenthetical comment, these are: 

1) The radio surface brightness is low by factor of 5-10 from other remnants of 
comparable diameter (This is fairly certain; Clark and Caswell, 1976). 

2) A low value of the blast energy E q is derived here. (Many assumptions used.) 

3) The supernova was not detected. This possibly suggests a low peak optical 
luminosity for the supernova. (Other plausible reasons exist.) Ttr*.s argues 
for a Type II supernova which is inconsistent with feature number 2. 

4) A central point source exists that is not a radio pulsar. (This is fairly 
certain.) 

How many of these are true anomalies and what possible connection exists 
between them is an open question. 

Data > 2.5 keV 

Figure 5 shows all the data displayed with the T e * T m model given in Table 3. 
The model underestimates the data at photon energies ^2,5 keV. We will restrict 
our discussion to the electron-ion equilibrium model since the "cold electron" 



model gives similar results in this energy range. The possibility exists that 
this anomaly is caused by an incorrect measurement of the background* The 
background of the SSS is knwon to have two components (Symkowiak, private communi- 
cation) : a "constant" background component that is simply proportional to the 

integration time of the observation and a "bump" background component that is 
also nominally proportional to the integration time but also has a scale factor 
relative to the constant background that can vary from one integration to another. 
I is customary to allow this scale factor to be a fitting parameter. 

Figure 6 shows the residuals at the end of the spectrum for observations 1, 

2, 3, and 4, respectively. The model used was the same one used to the fit the 
composite spectrum. The only variable in the input model between the different 
fits was the overall normalization factor which was determined only by the data 
with photon energies < 2.5 keV. The "bump background" scale factor was also a 
free parameter. Figure 5b displays the bump background component. The bump 
background alone cannot account for all the anomalous residuals. As can be 
seen, the shape of the anomalous residuals is different from the bump background 
component. In some cases the residuals have a maximum value at ^ 2.7 keV where 
the bump background is very small compared to its peak value. 

The inability of the standard background subtra^ f on procedure to account for 
the anomalous residuals does not preclude the possible 7 that the problem is 
related to background subtraction. Other variations in the intensity and spec- 
trum from the nominal background model are suspected in other SSS observations 
(Symkowiak, private communication) . Most of the background is produced by false 
counts made by charged particles in the magnetosphere of the earth. Variations 
in the density and spectrum of these charged particles can explain variations in 
the SSS background. Indeed, the "bump" background scheme was an empirical attempt 
to account for these variations but perhaps not comprehensive enough. 



Table 6 gives a measure of the residuals in each of the observations. The 
residuals were averaged between the energy of 2.5 and 3.2 keV. This range was 
chosen to minimize the dependence of measurement on the amount of "bump" back- 
ground. Also given, for comparison, is the average flux over the range of 0.8 
to 2.5 keV received during the various observations. The residuals from observa- 
tions 1 and 2 are twice as large as observations 3 and 4. This may be significant 
since the first two observations form a distinct set, grouped in time, from the 
last two observations. As can be seen in Table 1, observations 1 and 2 and 
observations 3 and 4 were close together in time but, the epochs for each pair 
were relatively far apart 9 days). The observations in the first pair looked 
at edges of the remnant and subtended a smaller fraction of ^he surface brightness 
than the second pair which was more centrally located on the remnant. The count 
rates from the first pair were smaller for this reason. The first pair also has 
shorter Integration times. The correlation of the intensity of the residuals with 
the epoch of the observation and the anti-correlation with the intensity of the 
known source intensity at other photon energies would support the notion that 
the residuals were caused by an Inaccurate model of the background. Anti- 
correlation of the residual intensity with integration time might also be indicative 
of background problems. Examination of other observations taken during this time 
have not been helpful. The source intensity at the pertinent energies overwhelmed 
the effect we were looking for in other observations. In any case, the effect 
might vary with position over the earth as well as the orientation of the telescope 
relative to the direction of the magnetic field and both of these parameters may 
be correlated with the pointing of the telescope in the direction of RCW 103. 
However, the evidence for background contamination is far from conclusive, and we 
yu examine possible contributions to these residuals that originate at the 


remnant . 



Attempting to explain the residuals in terms of emission from an optically 
thin plasma is difficult. The data require a rising photon spectrum at photon 
energies greater than 2.5 keV. Any additional continuum component, then, falls 
off faster than the data regardless of the temperature. Line emission is a 
possible explanation, but it has difficulties. Any significant contribution due 
to the line emission from elements with (Z _< 16) can be ruled out. The strength 
of the Ka emission features gives upper limits to the amount of power from these 
elements in the lines above 2.5 keV. For example, one of the most powerful lines 
of this type is a result ofann m 3ton=l transition in S XV (line energy at 
2.88 keV). Even under the most optimistic condition the power in the S Ka feature 
would constrain the power in the 2.88 keV line to be factor ^ 7 below that re- 
quired to explain the data at that energy. An overabundance of chlorine of a few 
hundred relative to solar could explain the data but is unlikely. 

Before considering other spectral forms, a relevant factor is the central 
point source of the remnant. Enough systematic uncertainty exists that it cannot 
be said for sure whether or not the point source was within the field of view of 
the detector during each observation. It was certainly within the field of view 
during observation number 4 and most likely within the field during observation 
number 3. Given the nominal positions, the point source was not involved in ob- 
servations number 1 and 2. Including possible offsets of the collimator (Symkowiak, 
private communication), it might be in one of those fields of view, but not both of 
them. It is not possible, then, for the point source to explain all the residuals. 
It is possible that the residuals in observations 1 and 2 are not caused by the same 
phenomenon as in latter observations, so we cannot rule out the role of the point 
source in this matter. 

A black body spectrum required a temperature of ^ 5 keV and an emitting region 
with a characteristic size of ^ 100 m at a distance of 4 kpc. Taking a neutron 
star density as an upper limit to the density of the emitting region, we calculate 
that the cooling time for this proposed black body would £ 100 days. Thus the 



emitting region would need to be ^ 0.01 the size of a neutron star and would need 
to be re-energized. 


Conclusions 

The main conclusions of this chapter are: 

1) No variation can be detected in the spectrum collected from different regions 
of the remnant. 

2) The data cannot be used to discriminate between whether non-Boulomb electron- 
ion energy excahnge processes may be present at the shock front. 

3) Assuming that non-Boulomb processes are present - a likely hypothesis given 
results from other young SNR's, the data are consistent with the idea that 
the emission is all from the shocked interstellar medium with approximately 
solar composition of heavy elements. 
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Table 3 

N1E Model without ejecta 



* 

j T, * T, 

best fit 

best fit 

t. = r t 

~ limits 

C, v (10 ,S cm -3 ) 

1.8 

056 

0.66 

0.35 

rj (foe l/s cm' 2 ) 

0.5 

0.34 

0.30 

0.49 

§ v a/q fo© 

• 

- 

- 

• 

r (fciloyears cm -1 foe -l/a ) 

3.6 

4.2 

4.8 

8.0 

Ng (10 21 cm' 2 ) 

5.0 

3.7 

4.1 

3.2 

Axt (solar) 

0.83 

1.5 

1.4 

2.1 

i4si (solar) 

0.79 

1.0 

1.0 

1.2 

Xs (solar) 

1.9 

1.3 

1.6 

1.3 

Are (solar) 

0.20 

0.75 

0.71 

1.3 

xi O' = 31 ) 

2.3 

1.9 

2.1 

2.1 

a (foe 1/3 cm) 

m 

0.37 

0.26 

01.1 

Eq (foe) 


0.05 

0.02 

0.6 

n 0 (cm -3 ) 


0.92 

1.2 

0.45 

f 3 (1C 9 years) 

0.85 

1.6 

1.2 

3.2 

D (-?c) I 

1.6 

2.6 

1.9 

6.7 


cb. foe s jo 31 ergs (fifty-one ergs; 0 Hans Bathe) 
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Table 4 

Two Temperature CIE Model 

Cxv.lrmi 

1.1 X I0 ,a cm* 5 

T 

*1 ,tow 

0.51 keV 

r - 

B.4 x 10 U cm* 5 

.high 

£6 keV 

Ny 

5.9 x lO 21 cm* 8 


2.3 solar 

A* 

1.9 solar 

As 

3.2 solar 

A Tm 

2.9 solar 
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Table 6 



Residual flux 
(2.5 - 3.2 keV) 


Total flux 
(0.8 - 2.5 keV) 


4.5 ± 0.9 
4.7 ± 0.9 
2.1 ± 0.6 
2.3 i 0.5 


26 5 ± 0.9 
28 6 ±1.0 
49.3 ± 0.8 
38.1 ± 0.7 


10" z cuts (SSS area)" 1 s" 1 keV" 1 









Figure Captions 


Figure 1 

The x-ray surface brightness map (Touhy, private communication) from the 
Imaging Proportional Counter (IPC) experiment on the HEAO A- 2 spacecraft. 

Circles represent the four fields view of the SSS observations. 

Figure 2 

The results of attempting to find a generalized model that was consistent with 
the four different observations are displayed. Histogram represents model. 

Figure 3 

Results are shown of best fits to the data when no ejecta was included (far Sedov 
limit) . (a) is for Coulomb heating model and (b) assumes electron-ion thermal 

equilibrium. Insert in (b) is the projection of X^ space on the r\ - x plane. 
Contours show the change of X^ from the best fit. Both 4.6 and 10 are shown. 

Figure 4 

Results are shown of best fits to the data when ejecta was included, (a) is for 
Coulomb heating model and (b) assumes electron-ion thermal equilibrium. 

Figure 5 

Figure 5 shows the same fit as was in Figure 3b but, nov, all the data is shown, 
(b) displays the residuals of the fit and (c) gives the bump fraction function 
for comparison. 


Figure 6 


The residuals from each observations to the model in Figure 3b are shown. 
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4 * 4 Cone lud ing Remarks 

The NIE models seemed successful from the point of view of the "global" 

issues since the models did give adequate fits to the data with physically 

2 

reasonable parameters. The x were somewhat high, but this is reasonable given 
systematic uncertainties in the modeling of the detectors, known uncertainties in 
the atomic physics, and deviations at the source from the ideal conditions 
assumed in the model. Although the existence of these systematic uncertainties 
gives reason for some caution in accepting inferred remnant parameters, the 
physically reasonable results do indicate that these models do hold some promise 
for use in future observations. 

CIE model fits to MSH 14-63 data had required two plasmas at different 
temperatures. The NIE model showed that with the exception of the S and Fe Ka 
lines, the data could be accounted for by shocked ISM of solar composition. The 
offending lines only required abundance changes by a factor of ^ 3, but as was 
discussed they could be explained by the density inhomogenities that do seem to 
be present in MSH 14-63 or by ejecta. 

In the case of the SSS data for RCW 103, CIE modeling again required two 
temperature components in order to give adequate fits. Abundance anomalies of 
£ 2 were reouired for Mg, Si, S, and Fe to explain the intensity of the line 
emission. The NIE once more proved successful in providing good fits to the data 
assuming emission from shocked ISM of near solar composition. 

The questioi ■ . emission from reverse-shocked ejecta was not clearly 
answered by the work here. For thc3e data sets, 6 was not a well determined 
parameter. Wnen 6 was low, the ejecta was too small an emission component. With 
increasing 8 the emission started to become significant but the temperature was 
so low that, the bulk of the flux from the ejecta component was either below the 
bandpass of the detector or was attenuated by interstellar absorption. 



On the issue of electron- ion thermal equilibrium the data from MSH 14-6' 
resoundingly asserted that some non-Coulomb energy exchange mechanism must exist 
behind the blast wave. MSH 14-63 joins Cas A and Tycho as young remnants where 
evidence for a non-Coulomb mechanism have been found. The case of RCW 103 is 
uncertain. Lack of information in the spectrum above 4 keV makes the interpreta- 
tion difficult. 

Future observations providing simultaneously good spectral and spatial 
resolution which are planned for the next generation x-ray telescope, can be 
invaluable for separating the ejecta component of the spectra from the shocked 
ISM and for establishing whether surface brightness variations within a remnant 
are due to a filling factors or density inhomogenities. Observations with spectral 
resolution of AE/E £ 1% with or without spatial resolution would allow measurement 
of lines from individual ionic species of the same element. These would allow 
more direct plasma diagnostics than those in this thesis which had to infer the 
fine structure from lower resolution spectra. 



5.0 Optical Identification of X-ray Sources 

Discovery of new X-ray sources is dramatic and can have considerable 
scientific impact, but until an optical counterpart for the X-ray source is 
found little further work can be done. This is true because X-ray detectors have 
comparatively low spatial and spectral resolution, and are encumbered by the 
necessity of space-borne observatories. Once an optical counterpart is found 
the full range of accumulated optical experience can be brought to bear, 
especially in the areas of high resolution spectroscopy and long and short-term 
variability. An optical position is also precise enough so that instruments in 
other wavelengths; ultraviolet, infra-red, and radio, for example; can be 
efficiently utilized. 

The specific set of X-ray sources we are concerned with are those discovered 

by the HEAO-i A-2 experiment Low Energy Detectors (LED’s). This experiment 

surveyed for the first time over 95% of the sky in the energy band 0.15-2.8 keV 

-11 -2 -1 

down to a limiting magnitude of 1 x 10 ergs cm s . A catalog of the sources 
seen by this experiment, using a very stringent 6a existence criteria, contained 
114 sources (Nugent et al 1983) . To appreciate the impact of this experiment 
consider that of the 49 sources in this catalog with flux between 0.15 and 0.44 keV 
only 12 were previously known. 

Unfortunately the uncertainty in position of these sources, resulting from 
the fluctuations in counting statistics was significant . The position was only 
determined to a characteristic size of ^ 0.2° in one direction and ^ 3° in the 
other. As the brightest objects in the X-ray sky are not necessarilly bright 
optically this is far to imprecise to allow optical identification. 

To further refine the position it is necessary to acquire more data. We have 
employed three techniques. In one, the Einstein X-ray Observatory was used to 
search the HEAO-1 LED error boxes for bright X-ray emitting objects. The Imaging 
Proportional Counter (IPC) aboard Einstein has a field of view of ^ 1° and 



produces X-ray source positions accurate to 1*. An object brighter than twelfth 
magnitude optically in such a small box was sufficiently uncommon as to make the 
identification by coincidence. Fainter objects had to be studied spectroscopically 
to select the X-ray emitter. 

The portion of this program partially supported by this grant was the optical 
identification. Utilizing the Mt. Palomar 60" telescope during the phase at 
Cal Tech and the Black Moshannon 60'* telescope during the phase at Penn State, we 
performed high and low resolution spectroscopy. 

Of eight HEAO-1 source error boxes completely searched by Einstein six 
produced strong X-ray sources, (The other two are presumably variable and so not 
detectable during the Einstein observation.) Two of these are located in the 
Southern Hemisphere and not accessible to lit. Palomar or Black Moshannon, but we 
have managed successful cooperation with our southern colleagues (Jensen, Nousek 
and Nugent 1982; Mason et al 1983; Nousek and Pravdo 1983). 

In addition to the bright HEAO-1 sources fainter X-ray sources appear 
serendipitously in the IPC fields. We have endeavore ' to identify these cs well. 

The observing pattern was as follows. For the bright stars, brighter than 
12 m , we assumed the star was the X-ray source and undertook high resolution 
spectroscopy searching for anomalies in the spectrum. For the fainter objects we 
conducted imaging color photometry to select blue candidates in the field and 
followed up with low resolution spectroscopy down to ^ 16 m . At Mt. Palomar we 
conducted 18 nights of observing, ten nights to spectroscopy, four nights to 
imaging photometry and four were clouded out. In this time we performed 
photometry in three colors on 18 X-ray fields, mostly seredipitous. Low 
resolution spectral types are obtained for 36 possible candidates, as well as ? 
modest library of comparison spectra. Repeated spectra of 18 br- candidates 
were obtained at high resolution to look for spectral changes with time. 



At Black Moshannon 31 nights were assigned, of wnich 18 were lost to clouds 
or equipment failures. The remaining 13 nights wore devoted to low resolution 
spectroscopy of 18 objects, one of which was observed 50 times to follow line 
profile shifts. 

The scientific return of this program has been substantial. A new Seyfe.rt 1 
galaxy was discovered at a redshift of Z = 0.038, E1615+061 (Pravdo et al 1981a, 
Pravdo et al 1981b). The presence of a Seyfert In a catalog of the brightest 
soft X-ray sources is surprising. The optical properties of E16154061 seem 
ordinary (for a Ceyfert). Pravdo et al suggest the great soft X-ray luminosity 
may be typical of Seyferts, but only rarely able to escape the core region of the 
galaxy. Hence only a few Se” ts at any one time will exhibit the X-ray properties 
of E1615+061. 

The apparently normal nintn magnitude G star HD155638 shows dramatic Ca H and 
K emission reversals (Nousek et al 1980, Stern et al 1981). These emission 
reversals are diagnostic of an active chromosphere, which apparently results in 
a powerful coronal X-ray source. The centroid of the reversal shifts with respect 
to the absorption profile but the radial velocity of the sharper lines suggest 
no orbital motion. Either HD155638 is a binary system of two extremely similar 
stars. In which only one is active, or it is a single star rotating rapidly with 
highly confined surface regions of activity. While HD155638 shows some similarity 
to the RS Canes Venaticorum stars, the most acti e coronal emitters, there are 
still some differences. 

A very recent discovery has been a new AM Herculis type binary system 
E2003+225 (Nousek et al 1982, Nousek et al 1983). These systems consist of a low 
mass star orbiting around a highly magnetized white dwarf. Rather than the 
accretion disk found in most astrophysical. accretion situations in AM Her stars 
the strong magnetic field channels the flow into an accretion columns. The 
study of these systems is particularly intriguing because of the complex behavior 



of the emission line profiles and doppler shifts, strong and variable optical 
polarization and high UV and X-ray luminosity and variability. Of about a dozen 
AM Her systems known E200 3+225 has the longest period (3 n 42 ) and thus both 
makes observations easier and probes the effect of system separation on the 
phenomena of the type. A paper combining HEAO-1, Einstein , high speed photometry 
and circular and linear polarization is currently in preparation. 

A similar systems with weaker magnetic fields, AC Cancri, is a strong 
candidate for another HEAO-1 source. Time resolved spectroscopic observation 
of AC Cnc through eclipse with the Palomar 200** telescope are currently being 
reduced . 

Useful as this technique is, the demise of the Einstein Observatory has 
forced us to look for alternative means of identifying the remaining 36 sources. 

Of the HEAO-1 sources located by Einstein nearly all were typified by emission 
line spectra, and all are relatively bright (£ 15 m ) . Optical techniques to 
search for emission line and blue continuum objects exist. One of the best is 
to use a wide field Schmidt ;.e]esrope with an objective prism. We have taken 
at series of 22 plates in th-Ls way , using the 24" Burrell Schmidt at Kitt Peak 
National Observatory and the 24" Curtiss Schmidt at the Cerro Tololo Inter- 
American Observatory. With a single exposure the Schmidt telescopes take low 
resolution spectra of every object within the 5° x 5° field of view. This is 
easily large enough to include the entire HEAO-1 error box. The major problem 
lies in the analysis of the plates. At galactic latitudes below 20° the confusion 
of spectra prohibits the use of this technique. 

The plates token in these two runs (April and November, 1982) have been 
searched and seven promising candiates resulted. Higher resolution spectroscopic 
work and IUE observations are underway to better identify these objects. 

The third optical technique we are using to study X-ray souces is to search 
for rapid variability. Jensen (1982) has used the HEAO-1 LED data to contrain 



the X-ray emitting properties of the old nova and mm like variable sub-class of 
cataclysmic variable. Current theory suggests these variables are only in out- 
burst infrequently, but the observed number of outbursts is relatively high. If 
this is true then many of these variables arf near us 'waiting* to go into 
outburst. Their X-ray emission should be statistically detected as a population 
of moderately bright X-ray source. 

Jensen (19*2) finds there are too few such source; to resolve the missing 
dwarf nova problem and concludes the X-ray emission must oe smothered. If this 
is true then the old novae can only be detected by their rapid optical variability. 

In order to search for variable stars, shore, 15 m ) exposures using the 

Palomar 48" Schmidt telescope have been taken of the same 7° x 7° part of the sky. 

A series of nine such plates, containing 27 images, have been collected. Most 
stars will produce 27 approximately equal images, but stars which vary over 
short timescales will produce varying images. 

The major difficulty in this program is selecting out the varying images 
from the v 10** star images on each plate. The Image Processing Lab at the 
Jet Propulsion Lab has kindly offeree to digitize and analyse a trial plate. 

This cc ;>uter selection process is not oerfect, but from the top 50 candidates 
picked by the computer we have found evidence of variability in six on tl esc«les 
of 30 minutes or less. Follow up spectroscopic studies will be undertaken at 
Black Moshannon Observatory this summer. 



6.0 Evaluation of certain Metal Oxides as X-ray Detectors 

Current X-ray detectors possess energy resolution which is determined by the 
fluctuation in the number of ion-electron pairs generated in the detector 
following photoelectric absorption of an X-ray. For gaseous detectors the energy 
required to produce an ion-electron pair is of the order of 25 eV. The fluctuations 
in the signal generated by a mean number of ion-electron pairs is proportional to 
tb. square rot of this number for a given material times a factor (Fano factor) 
wl \ch modifies this statistical process based on the production of ion-electron 
pairs by the photo election. Counter effects, such as variable gas gain and 
charge collection, are sometimes lumped with this correction to provide an overall 
correction to the statistical fluctuations in the mean charge generated in the 
interaction. 

Solid state detectors based on Si, Ge or Hgl^ require less energy to generate an 
ion-electron pair than a gas, since the valence and conduction electrons require less 
energy to remove them from the nuclear charge when the atoms are bound in a solid. 

The energy for these solids are roughly 4 eV to generate an ion-electron pair. 

Transition metal oxides are thought to require less energy to produce an 
ion-electron pair than other solids because of the nature of d-electrons forming 
the bonds in the crystal lattice. The energies may be as low as one eV or 
slightly more. If such a material can be found for which electron trapping is 
not high, then a factor of up to two in intrinsic energy resolution could be 
achieved. No direct measurements of this quantity has been made on transition 
metal oxides, so a program was initiated to obtain some basic data on several such 
materials including ^2°s 9 Cr 2°3* Ti0 2 and Mn0, of which stlow e l e ^ tri cal 
properties at low temperature which can yield data on the ion-electron pair pro- 
duction and electron trapping. 

A crystal of each of the above substances has been obtained. A laboratory set 
up has been developed for coating the crystals with gold to provide electrical 



contact. A dewar for cooling and a thermistor thermometer has been constructed 
for the measurements. Unfortunately, the graduate student working on this 
problem decided to go into theory, so the measurements have not been completed 
at this time. They will be carried out by a lab assistant In the next few 
months. 
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